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Riassunto. — Si applica al calcolo dell’energia dello stato fondamentale 
dell’ H? e dell’Het un metodo assai più rapido e preciso dei precedenti. 
Esso & basato sull’idea di trovare il minimo del valor medio dell’ Hamil- 
toniano, non tra una classe di funzioni contenenti un certo numero di 
parametri, ma tra l’intera classe delle funzioni di un determinato argo- 
mento (che viene assunto essere la somma dei quadrati delle distanze 
tra i nucleoni). I risultati e le conclusioni sono esposti nel n. 5. 


La presente nota & un contributo al problema dell’energia di legame del- 
1’H: e dell’Het; in essa viene descritto un nuovo metodo di calcolo dell’energia 
di legame di tali nuclei, metodo che presenta sopra i precedenti vari vantaggi. 

Tl metodo in questione pud essere senz’altro applicato al caso in cui le 


ue nucleoni siano non centrali; tuttavia per semplicita in 


forze agenti tra d 
non ¢’é bisogno di 


questo primo lavoro considereremo soltanto forze centrali: 
insistere sul fatto che le forze tra due nucleoni contengono in effetti una grossa 
percentuale di forze non centrali le quali contribuiscono, per lo meno nella 
forma attuale, considerevolmente ai valori dell’energia di legame (1). I cal- 
coli che presentiamo hanno percid più che altro lo scopo di illustrare gli accen- 


nati vantaggi del metodo. 


(1) E. Gerrouy e J. Scuwinger: Phys. Rev., 61, 138 (1942); Tst-Mine Hsu e 
H. S. W. Massey: Proc. Roy. Soc. (A) 196, 135 (1949); R. E. CLapp: Phys. Rev., 76, 
873 (1949); Tsı-Mıng Hsu e Kune-Neou Hsu: Phys. Rev., 78, 633 (1950). 
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Fissando l’attenzione sulla espressione della energia potenziale tra due 
nucleoni: 


a 


(1) Hira) = — At (2) {oo + OP, + APP, + MP.) : 


assumeremo per f una gaussiana: 
2 
exp [— rx/a?] . 


Il metodo puö tuttavia essere egualmente bene applicato ad un potenziale yuka- 
wiano od esponenziale. E fisseremo l’attenzione sulla seguente scelta dei para- 
metri: 


A = 35,6 MeV a — 2,25-10- cm 


g =b+h= 0,19 (m+w+b+h=1). 


La ora detta scelta della forma dell’energia potenziale e dei parametri & 
stata fatta per poter paragonare i nostri risultati con quelli di vari altri 
autori (?). 

Nel n. 1 accenneremo ai risultati dei vari lavori finora pubblicati con la 
scelta (2) dei parametri; nel n. 2 esporremo il metodo per l’H*; nel n. 3 per la 
particella «; nel n. 4 raffineremo il calcolo dell’H? e infine nel n. 5 discuteremo 
i risultati ottenuti. 


1. — Situazione del problema. 


La situazione del problema é illustrata dalla tabella seguente: nella prima 
colonna di essa si é riportato il procedimento di calcolo usato, nella seconda 
la funzione d’onda di prima approssimazione usata, nella terza il numero delle 
approssimazioni successive (il termine « approssimazioni successive » avendo 
naturalmente significato diverso in un metodo variazionale, variazionale pertur- 
bativo, iterativo, ecc.), nella quarta il risultato globale ottenuto prescindendo 


dall’energia coulombiana, come faremo in tutto il seguito. 


(2) H. MARGENAU € D. T. Warren: Phys. Rev., 52, 790, 1027 (1937); H. MARGENAU 
e W. A. TYRRELL (Jr.): Phys. Rev., 54, 422 (1938); N. SVARTHOLM: The binding energies 
of the lightest atomic nuclei (Lund, 1945). 
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Funzione Numero | 
Metodo d’onda di ap- | approssi- Risultato Ref. bibliografica | 
| usato prossimazione | mazioni | globale | e annotazioni | 
| . | | 
zero successive) —E MeV 
variazionale gaussiana 13 7,21 MARGENAU e WARREN (?) | 
| 
: : ; £ | 
HS iterativo gaussiana 2 7,3 SVARTHOLM (?) 
F | LE fe 
presente lavoro — 2 7,79 (senza il raffinamento di | 
cui al n. 4 si ha — 7,24 
MeV) 
en | —_— eS a 
variazionale gaussiana 14 26,67 MARGENAU e WARREN (?) | 
(coordinate) 
individuali) 
variazionale gaussiana 11 26,76 MARGENAU € TYRRELL (?) 
Het (coordinate 
relative) 
iterativo gaussiana 2 ~27,1 SVARTHOLM (?) | 
presente lavoro _ 1 26,77 —- 
| 
| 


Nell’ultima riga della tabella si trovano i risultati del nostro lavoro. 


2. — Impostazione del metodo e applicazione all’H°. 


Tl metodo che useremo si ritrova già in un vecchio lavoro di HYLLERAAS (°) 
sull’atomo di He. Con tale metodo HyLLERAAS giustificava il fatto che una 
buona approssimazione ad una corretta funzione d’onda per lo stato fonda- 
mentale del parHe si poteva gia ottenere con una funzione del tipo: 


Fe 
exp = 


0 


(+ rl ) 


rappresentante gli elettroni muoventisi in un campo coulombiano schermato. 
Noi abbiamo semplicemente applicato tale metodo ai problemi nucleari 
come ora più in dettaglio passiamo a vedere. 
Fisseremo l’attenzione su una funzione d’onda del tritone antisimmetrica 


(3) E. A. Hytupraas: Zeits. f. Phys., 54, 347 (1929); 48, 469(1928). 
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di spin rispetto allo scambio dei due neutroni e completamente simmetrica 
spazialmente rispetto allo scambio dei tre nucleoni; ossia scriveremo la ® 
dello stato fondamentale dell’H® sotto la forma: 


(3) C= er % (G33 — %ß,)p(xı , Xe, 3) ; 


dove «x, indica la funzione di spin del protone, &, Ps, &, ß, quelle dei neutroni, 
%1, %, x; le coordinate dei tre nucleoni e la p é simmetrica rispetto à x,, x2, x3 
e normalizzata. E noto come una funzione d’onda della forma (3) avrebbe la 
simmetria corretta rispetto a permutazioni se non vi fossero nella energia 
potenziale i termini di Bartlett e di Heisenberg (prima approssimazione di 
Wigner). In effetti con la presenza di tali termini, la funzione d’onda corretta 
dello stato fondamentale viene ad essere una combinazione lineare di due fun- 
zioni, una delle quali ha una parte spaziale simmetrica rispetto allo scambio 
di x, con x,, mentre l’altra ha la parte spaziale antisimmetrica. Nell’usare 
perciö come funzione variazionale una funzione del tipo (3) commettiamo 
certamente un lieve errore; tale errore & perd notoriamente del tutto trascu- 
rabile e non ce ne occuperemo più, essendo gia stato valutato da vari autori (4). 

Introducendo allora coordinate normali relative e baricentrali q, r, X 
legate alle x,, x, x; da: 


x x 3 il 
(4) an, ee (x; — xs), Pa (ar a a) 


e, sottraendo l’energia del baricentro, il valor medio dell’energia dell’H* nello 
stato (3) si scrive: 


2 


| exp 


„2 


— =) p dt, dt, , 
a? 


Fa [ es 
1e 7 AM al A) exp 


neutroni 


dove @ € 7ix, Ya; Si intendono ora espressi mediante le r, q e si è posto: 
(5) ‚Br (1-35) A, C — (1—29)A. 


In uno stato S come quello con cui abbiamo a che fare la funzione d’onda œ 
del sistema dipenderà semplicemente dalla configurazione del triangolo costi- 
tuito dai nucleoni, configurazione che pud essere caratterizzata dal modulo r 
di r, dal modulo q di q e dall’angolo y tra reg. 


(4) W. RARITA e R. PRESENT: Phys. Rev., 51, 788 (1937); H. MARGENAU e W. 
TYRRELL: Ref. (2). 
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Posto 
a 4M 4M 4M 


B= By 0 E* 


il minimo valore di E* &, per uno stato 8, dato dal principio variazionale : 


: : dp\?. [àp\ dp\? 
© aff [an {an re (Be een 


15 
Saree 


3 2 
— sin y gr? (Bs > exp - a = O* exp 
2k 


con la condizione 
CO co va 
(7) far | aa | sin y dy rate? il 
0 0 0 


E a questo punto che il nostro procedimento si distacea da quelli conven- 
zionali per riattaccarsi all’idea usata da HYLLERAAS per Vatomo d’He. 

Noi cercheremo precisamente una funzione variazionale che soddisfi alla 
(6) (7) tra la classe delle funzioni del solo argomento completamente sim- 
metrico 9? = r?- q’. 

Si vede immediatamente che &: 


(8) 1 = 2 SE 2% are xs — (%1°X2) — (= (X %3) » 


e l’espressione a secondo membro della (8) non é altro che la meta della somma 
dei quadrati costruiti sui lati del triangolo eostituito dai tre nucleoni. Intro- 
ducendo coordinate o, 0: 


(9) r= ocosé, q=esind, 


il principio variazionale (6) (7) si scrive allora: 


co 7/2 
(10) | ao | dd +205 cos? 0 sen? 0: 
0 0 


Oe N ee le 
ee (2B* + C*) exp |— 0 al OF 


(avendo gia compiuto l’integrazione rispetto a y tra 0 e x) con la condizione 


x/2 


co 
(11) 2 | de | aoe cos? 0 sen? 6g? = 1. 
0 


0 
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Nelle (10) (11) la @ é intesa ora come funzione della sola o; la notevole sem- 
plificazione rispetto alla (6) à dovuta al fatto che la y & simmetrica, cid che 
consente di scrivere nell’energia potenziale con simboli evidenti: 


(12) [asp = [ (13) 9 =| (23) +. 

Gli integrali rispetto a 9 nelle (10) (11) possono essere fatti e l’equazione tipo 
Schrodinger per la p(o), equivalente al principio variazionale in una sola varia- 
bile al quale le (10) (11) danno luogo quando si compiano le dette integra- 
zioni, viene ad essere: 


d°p YP é 4Ma 0? : 202 2M 
(13) apna hc ae he 2 2! (ie 2 Ep =0, 
con la condizione che 
CO 
(14) [or do, 
0 


sia finito il o che compare nelle (13) (14) è pari a 4/2/4/3 volte il o che con- 
pare nella (10)). 

La determinazione del più basso autovalore della (13) fornisce dunque per 
quanto abbiamo detto, il migliore valore dell’energia di legame dell’H® che si 
puö ottenere con una funzione della sola o. Nella (13) J, & la funzione di Bessel 
di ordine 1, quale definita ad esempio in (5). 

Posto 


(14) oe’ ®plo) = F(p), 


e introducendo la variabile adimensionale: 


(15) er 
a 
posto inoltre: 
2Ma? 
(16) E£ — ar -E y 
la (13) si serive 
4 Ma? exp [— x?] 15 


17) F’+/|(2B+ 0) { — iJ, (iw?) }— 


je a? 4x? 


(5) W. MaGxus e F. OBERHETTINGER: Formeln und Sätze für die speziellen Funktionen 
der mathematischen Physik (Springer Verlag, 1948). 
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Il grafico di quella che possiamo chiamare l’energia potenziale efficace 
penaläEts: 


15 4Ma? exp [— 2?] 


re ar = {— ii (ia®)} , 


(18) V(x) = 


è riportato nella fig. 1 (curva continua). Si tenga presente che coni valori (2) é: 


h? 4Ma? 
Maya 8,1 MeV, (2B + 0) —— = 42,4548 . 


(9) h? 


La salita della curva nell’origine & da mettersi ovviamente in relazione con 
il principio di indeterminazione. 
Risulta dalla figura come V(x) abbia la forma di una energia potenziale 


H3 


V(x) 
ee I Een 1010) Fo (x) 


\ 
1 

\ 

\ 
CE tt tt tt 
01 02 03 04 | 06 07 0809 | I) 12 13 1415 LS 17 18 19 2 22 23 


El en ee te 


per molecole. À facile vedere che essa pud essere assai bene approssimata da 
una energia potenziale alla Pöschl-Teller introdotta (*) da questi Autori per 
studiare le vibrazioni delle molecole biatomiche. Tale energia potenziale, i cui 
autovalori sono ottenibili esattamente, ha la forma generale 


il A B 


à we, Ge 
(20) 12 (2) a2 |sinh? « cosh? «x 


A,D,a 


in cui A, 2, « sono tre parametri arbitrariamente fissabili. 


(6) G. Pöscut e E. TELLER: Zeits. f. Phys., 83, 143 (1933); vedasi anche: L. INFELD 
e T. E. Hutz: Rev. Mod. Phys., 23, 34 (1951). 
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Scegliendo: 
(21) HE, DEI, En 


si ottiene la curva tratteggiata in fig. 1; come si vede tale curva approssima 
assai bene la V(x) ed anche nell’origine la differenza tra (18) e (20) è finita 
(ed assai piccola — 2,27 MeV). 

Gli autovalori e del potenziale di Pöschl-Teller, posto: 


(22) u(u+1)=S8, u> 0 
vy—1)=Z, y > 0 

sono 

(23) EP) = — au — y — 2n)?, 

con: 

(24) ea ee 


À 


e l’autofunzione normalizzata dello stato fondamentale (n = 0) é: 


_ 2al(u+1/2) |" 
Io + 1/2)T(u — ») 


(25) io. = sinh” «x cosh~" ox. 


L’autofunzione (25) relativa alla scelta (21) delle costanti verrà indicata nel 
seguito con (x). Con la scelta (21) di H, 2 ed «, si ottiene » = 2,5000 
e u = 3,8874 e quindi: 


3 h? 3 > mm 
(26) CEE DE 1, MeV, 


L’autovalore vero & della (17) (che d’ora in avanti affiggeremo di un indice (#9) 
si pud ottenere con ottima approssimazione (e in ogni caso con approssimazione 
sicuramente per eccesso), considerando la differenza tra V(#) e P(x) come per- 
turbazione; l’integrale: 


co 


(27) so) { V(x) — Pisya, 10, 1 (a) } da, 


0 


é stato eseguito numericamente e vale: + 0,136 di modo che & certamente 


he 
(28) ŒE = Ma (H) ¢ << — 7,79 + 4,05-0,136 = — 7,24 (MeV). 
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Dalla disuguaglianza (24) risulta inoltre che nell’H® si è lontani dall’avere 
stati eceitati S stabili. 


3. — Applicazione alla particella «. 


Lo stesso metodo usato per l’H* pud essere trasportato alla particella «. 
Anche qui la parte spaziale della funzione d’onda & con ottima approssima- 
zione completamente simmetrica. Useremo in questo caso un procedimento 
alquanto abbreviato, supponendo sin dal principio che la funzione d’onda spa” 
ziale dipenda soltanto dalla somma R? dei quadrati costruiti sui lati del tetra- 
edro costituito dai quattro nucleoni. 

Introdotte coordinate relative e baricentrali 


11/3 
ee T2 + Recs la) > 
si pud controllare immediatamente che è: 
1 2 2 2 
(30) Gee eae Ua aie Os 


Posto ancora: 

(31) Qie=t,0080,, Qw = t sind, COS; Qi. = 4 Sin 0, Sing, , 

e analogamente per Q, e Q,, si ha: 

(32) R=VE+R +; 

limitandoci a funzioni variazionali @ della sola R conviene ancora scrivere: 
(33) t, = Rsené cos 7 , Dr OS tj = Rsenésiny, 


(con il che 0<R< 00, 0 —E<x/2, 0 << x/2); ed il principio variazionale 


470 G. MORPURGO 


si serive allora: 


CO 7/2 a /2 
(34) 6 le an | sine € cost Ea [cost sin® yay 
0 0 0 
dp \? _ = 4R? cos? 
= ll = 4 2 I) An) : = 
{ (24) — 4B + 20)- exp FEES) 0, 
dove: 
— 4M a UIE = 4M 
B SSS } — — y E == . 
Hani © h En eo 
con la condizione: 
co la ul? 
(35) [rar | sind € cos? Edé | cos? 7 sin? 7 An = 1. 
0 0 0 


Anche qui le integrazioni rispetto a € e 7 sono effettuabili e, come gia nel caso 
dell’H’, si perviene ad un principio variazionale nella sola variabile R (che 
non stiamo a scrivere) il quale risulta facilmente essere equivalente alla seguente 
equazione: 


d?p 8 do 4M u R 4M 

ae ER 20 = ee) 
er dr Rak? are 
dove: 

R\ 105 (I(k) 15 3 au 2R 
reg el 1 | X = 2 i 
@ (7) 166 |e | + dee aa] OPE PV ope ae J in: 
e 
iy 
ZA) = [exp [— V2]aV, 
6 
eS 
con la condizione che [pr dk sia finito. 
0 
Posto 

(38) ok = F, 


e introdotta oltre a cid la variabile adimensionale 


(39) EN. 
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posto ancora: 


2Ma? 
(40) eg à 
pe 

si ha l’equazione: 

2Ma®° | æ 12 
(41) Fr+ | (4B +20)8 ()- = F4+ eF=0, 

he A) m 
co 


con la condizione che | F:dx sia finito. L’energia potenziale efficace in questo 
0 


caso è: 


eee eee SL ate = 
S(@) == Re (4B + 20) @G ray 


ed & rappresentata (curva continua) nella fig. 2 in funzione di x. Anche qui il 
potenziale di Pöschl-Teller (20) con la seguente scelta delle costanti: 


5 1 
(42) Z = 12, 3=6l,5, HE FEU 


Ha } Het 


0 LS FE 
O1 92 03 04 05 


Fig. 2. 


approssima bene (vedasi la curva tratteggiata in fig. 2) l’energia potenziale 
S(æ) e poiché ai valori (42) corrisponde > = 4,000, u = 7,358 risulta 


Be 
2Ma? 


(43) (EP) = (a) e(P) == 97,022 MeV, 


472 G. MORPURGO 


la funzione d’onda corrispondente è: 


% x 
(a) 7) ~ Si 4__ eogh-7-358 —_ 
Fox) ~ sinh 13 osh 13? 
L’integrale: 
ce 
(45) | (‘OF o(x))” { S(æ) Jp 61,5, 1/1,3 (a) } da, 
0 
calcolato numericamente vale: + 0,105, onde risulta: 
84.05 
(46) () BE < — 27,022 + 0,105:—— = 26,77 MeV. 
1,69 
Contrariamente a quanto avviene per l’H* l’energia potenziale S(x) con- 
sente qui la presenza di un livello eccitato 2S con energia totale negativa. 
Il fatto che esistesse lo stato eccitato S era gid stato notato da MARGENAU 
e TYRRELL (?) i quali avevano trovato per l’energia di tale livello —1,3 MeV. 
La (23) e la (42) danno per il livello in questione: 


he 1 
2 Ma? 1,69 


(47) (4—v— 2)? — — 4,41 MeV. 


\ 


Naturalmente questo valore non è corretto; l’energia di questo stato eccitato 
e in effetti maggiore del valore (47) perchè l’autofunzione °F, del potenziale 
di Pöschl-Teller corrispondente al livello eccitato in questione non & ovvia- 
mente ortogonale all’autofunzione dello stato fondamentale di S(æ). Non insi- 
steremo su questo punto (sebbene sembri abbastanza chiaro che un valore 
assai migliore di — 1,3 MeV sia ottenibile (7) per l’energia del livello 28), perchè 
lo stato in questione non presenta alcun interesse, in quanto instabile per 
disintegrazione in due deutoni oppure in un tritone (He?) e in un protone 
(neutrone). 


4. — Raffinamento dell’approssimazione nel caso dell’H?. 


I valori dell’energia dell’H? e dell’He* che abbiamo ottenuto nei numeri 
precedenti basandoci sull’idea di cercare il più basso autovalore associato alla 


(7) Per fare un calcolo accurato (e in ogni caso errato per eccesso) dell’energia di 
legame de! considerato livello 28 basterebbe prendere come funzione variazionale 
CF, + 037, e considerare la radice più alta dell’equazione che si ottiene eguagliando 
a zero il determinante ottenuto minimizzando il valore medio deli’hamiltoniana rispetto 
a 0, C,. 
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classe di funzioni della sola 9 0 R possono essere naturalmente migliorati intro- 
ducendo nelle funzioni variazionali le altre variabili. 

Tale estensione sara qui fatta nel solo caso dell’H*; si ha di qui un’idea 
dei risultati che si possono ottenere per l’Hef. 

Una funzione completamente simmetrica nei tre vettori x,, x, x3, la quale 
sia, quando espressa mediante r, q, X, indipendente dal vettore X (e sia inva- 
riante rispetto a rotazioni), pud sempre essere espressa come funzione di tre 
argomenti completamente simmetrici indipendenti. Tali argomenti sono rica- 
vabili in base alle stesse caratteristiche geometriche del triangolo costituito 
dai tre nucleoni, e sono, ad esempio, oltre alla grandezza: ‚0? = 72 + q* pre 
cedentemente introdotta, l’area A del triangolo e il prodotto p delle lunghezze 
dei tre lati del triangolo stesso. Indicando con : proporzionalita, e riportandoci 
alle notazioni del numero 2 é: 


A : rqsiny : sin26siny, 
il 2 
DOTE (a 3 ) > cos? 6(3 — 4 cos? 0)? : cos? 36; 


la funzione d’onda più generale sara percid una (0, sin 29 sin y, | cos 30|). 

Per vedere l’influenza delle nuove variabili, faremo i calcoli con la funzione 

variazionale : 

(49) plo). + B cos? 30]: 

inserendo la (49) nel prineipio variazionale (6) (7), con calcoli simili ai prece- 

denti, seppure un po’ pit complicati, si giunge per la F(a) (che é collegata alla 

g(e) dalle (14) (15)) all’equazione (50) nella quale compare ora anche il para- 

metro ß che pud essere variato cosi da render minimo il più basso autovalore 

dell’equazione (50) stessa: 

2Ma? 
hi? 


(50) F"— V*(x)F + EF=0, 


con 


9 p? 1 42,4543 exp — a}. 
2 0.02 Oo = 


5. ; In. 2 
— tJ (iv?) — 27J (tx) + 5 xd (ie?) + = iyJ (1x*)| , 
ove Jo, Jay Js, J, sono funzioni di Bessel e le costanti o, t, y sono date da: 


3 
let, 


2 ß? 
De Ran 
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È chiaro che ponendo B = 0 nella (50) si ricade nella (17). À facile rendersi 
conto calcolando V*(#) per vari valori di © scelti nella regione in cui l’auto- 
funzione (I) (v7) à massima che la quantita V*(x) & in tale regione general- 
mente minima per un valore di # intorno a 0,22. L’integrale: 


co 


[gay { Pla) — V*a)} dr, 


0 


è stato valutato per f — 0,22 e si & visto avere un piccolissimo valore positivo 
dell’ordine di 10°, 

Segue che il più basso autovalore del potenziale V*(x) coincide entro le 
prime due cifre decimali con l’autovalore pitt basso del puro potenziale di 
Pöschl-Teller e vale quindi — 7,79 MeV. Di modo che l’energia dello stato 
fondamentale del’H* & certamente minore di — 7,79 MeV. 


5. — Discussione. 


Il metodo esposto conduce, come si sara osservato, a degli ottimi valori 
dell’energia di legame. Vogliamo qui elencare i principali risultati del metodo. 


1) Esso fornisce, all’approssimazione zero, un valore dell’energia di le- 
game altrettanto buono di quello ottenuto variazionalmente da MARGENAU 
e TYRREL e da MARGENAU e WARREN per l’H? e l’He? con un sistema di rispet- 
tivamente 14 e 11 funzioni; e, grosso modo, altrettanto buono di quello otte- 
nuto da SVARTHOLM; inoltre esso permette di valutare immediatamente, 
senza ulteriori calcoli, la presenza di livelli eccitati, aventi la stessa simmetria 
dello stato fondamentale. 


2) Il raffinamento (49) della funzione d’onda per V’H? conduce a un 
valore dell’energia per questo nucleo assai pit basso di quelli finora ottenuti; 
esso à tale da far ritornare in discussione la questione della coerenza dei para- 
metri sollevata da RARITA e PRESENT (4) (8), questione che tuttavia non. discu- 
teremo in quanto: a) sarebbe necessario svolgere gli analoghi calcoli raffinati 
per ’He?; b) & un po’ sorpassata per la indubbia presenza delle forze tensoriali. 


3) Conduce alla forma asintotica corretta delle funzioni d’onda, cid che 
è importante per trattare problemi di fotodisintegrazione e di fotoproduzione 


di mesoni; e la forma delle funzioni di onda & relativamente semplice (ve- 
dasi (25)). 


(?) L. ROSENFELD: Nuclear Forces (Amsterdam, 1948), pag. 305. 
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4) Giustifica parzialmente il «two equivalent body method» di FERN- 
BERG ($) (°) in quanto riconduce la risoluzione del problema dell’H? e dell’He* 
alla risoluzione di una equazione del tipo del deutone. 

Notiamo ancora che i calcoli numerici sono estremamente semplici e, per 
quanto riguarda Vapprossimazione zero, eseguibili in una decina di giorni. La 
maggiore semplicità dipende dal fatto che, una volta trovata una curva che 
approssimi bene il potenziale efficace, si tratta soltanto di calcolare integrali 
del tipo (27) (45) e che, sebbene nel presente lavoro essi siano stati calcolati 
esatti a tre cifre decimali, & possibile per tali integrali una rapida valutazione 
più grossolana, in quanto essi portanto soltanto una correzione relativamente 
piccola all’energia totale. Non si ha a che fare come nei metodi variazionali 
usuali con una situazione in cui l’energia totale si ottiene come piccola diffe- 
renza tra un grosso contributo dovuto alla energia cinetica ed uno cirea altret- 
tanto grosso dell’energia potenziale, di modo che il piccolo errore relativo nella 
valutazione di una di queste grandezze porta un forte errore sul valore ottenuto 
per Venergia. 

Vogliamo infine notare che in un recentissimo lavoro (gentilmente segna- 
latoci dal prof. MÖLLER) IRWING (°°) ha calcolato l’energia di legame dell’H? 
e dell’Het con un potenziale di Yukawa ed esponenziale facendo uso di fun- 
zioni d’onda che, nelle nostre notazioni hanno la forma: 


orexp[-ael e Rrexp —e«R], 


(rispettivamente per l’H° e l’Het) e minimizzando rispetto a x e ad n. 

Le funzioni variazionali usate da IRWING, che sono notevoli in quanto tra 
tutte quelle precedentemente usate in problemi nucleari sono le prime a godere 
di una certa semplieitä, pur avendo la corretta forma asintotica, rientrano 
nella classe g(o) 0 p(R) da noi usate. ; 

L’applicazione del presente metodo a nuclei superiori alla particella « sembra 
difficilmente possibile; il problema di cercare nel caso di un nucleo superiore 
alla particella « la migliore funzione tra la classe delle funzioni di un argo- 
mento convenientemente scelto si riconduce in tal caso, poiche a causa del 
principio di Pauli Vargomento & necessariamente non simmetrico, anziche alla 
risoluzione di una equazione differenziale ordinaria alla risoluzione di una com- 
plicata equazione integro-differenziale (1). 


Ringrazio il prof. B. FERRETTI per le varie utili discussioni avute. 
Desidero inoltre ringraziare, nelle persone del suo direttore prof. M. PICONE 


(°) H. Berue e R. BACHER: Rev. Mod. Phys., 8, 82 (1936). 
(10) J. IRWING: Phil. Mag., 42, 338 (1951). 
(4) Vedi ad es. J. A. WHEELER: Phys. Rev., 52, 1083 (1937). 
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e dei suoi collaboratori prof. A. GHIZZETTI e dott. W. Gross, l’Istituto per le 
Applicazioni del Calcolo per l’esecuzione dei calcoli numerici relativi alla valu- 
tazioue degli integrali: (27), (45), (52). 


APPENDICE 


a) Le formule usate nel testo per il calcolo dei vari integrali sono: 


n 
a2 


N DT ae! SLEOL ANS 
[sm 6 cos 00-58 2 ’ 9 i. Bim, n) — T(m + n) ’ 
0 

inn 


IAG) = i exp [iz cos t] cos nt dt ; 


2n 


de) EI eae) as J,(2). 


b) Si puö fare un controllo diretto della (10) inserendo in essa al posto 
di (0), exp [— ke?]. Tenendo presente la formula: 


ie a, (2b I + 1/2) 
| exp L- «419, (bi at = en 


0 


(Rey >0; Rea>|Imb)), 


si ottiene naturalmente il risultato (118) di BETHE e BACHER. 


SUMMARY 


A method much more effective than the former ones has been used to calculate 
the binding energies of H? and Het. The idea on which the method is based is to look 
for the minimum of the mean value of the Hamiltonian not among a class of functions 
containing some parameters, but among the whole class of the functions of an argu- 
ment (which is assumed to be the sum of the squares of the distances between the 
nucleons). With a gaussian potential and central forces (with the usual choice of para- 
meters of Margenau) one obtains in the first approximation — 7.24 MeV for H® and 
— 26.77 MeV for Het. These values have to be compared with — 7,21 MeV, and 
respectively — 26,76 MeV obtained by Margenau making use of a system of 14 an 
11 variational functions. A refinement for H* lowers its energy to — 7.79 MeV. Other 
advantages of the method (among which we mention the right asymptotic form of 
the wave functions obtained) are listed in section 5. 
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Sui metodi di soluzione 
delle equazioni della cascata elettrofotonica. 


F. MONTICELLI 


Assieurazioni Generali - Venezia 


(ricevuto il 30 Marzo 1952) 


Riassunto. — Le equazioni fondamentali della cascata elettrofotonica 
vengono ricondotte ad una equazione integrale di seconda specie del 
tipo di Volterra. Si dimostra che le soluzioni di SNYDER e di Scorr e quella 
di BABA e CHAKRABARTY soddisfano a tale equazione. 


Recentemente H. SNYDER (1) e W. T. SCOTT (?), con metodi equivalenti, 
il primo usando antitrasformate di Mellin, il secondo usando le trasformate 
di Mellin e di Laplace, hanno dato una medesima soluzione alle equazioni 
fondamentali della cascata unidimensionale elettrofotonica nell’approssima- 
zione B di Rosst e GREISEN (*) con sezioni d’urto asintotiche. ; 

Tl confronto numerico dello spettro integrale di questa soluzione con quella 
ottenuta precedentemente da H. J. BHABHA e S. K. CHAKRABARTY (*) ha rive- 
lato delle differenze che non sono a priori giustificabili (°). 

Scopo di questo lavoro é di ricostruire tali soluzioni per altra via, partendo 
cioè da una equazione integrale di Volterra di seconda specie, ottenuta con 
opportune integrazioni delle equazioni fondamentali nello schema della trasfor- 
mata di Laplace. 

Questo attacco si à dimostrato assai fecondo in quanto ci ha dato imme- 


31 - Il Nuovo Cimento. 
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diatamente l’equazione di partenza di SCOTT e successivamente la sua eq. (10). 
Inoltre, ci ha consentito di usare una rappresentazione in serie di potenze 
di (— ß/E), che ci ha condotto immediatamente alla soluzione (14) di BHABHA 
e CHAKRABARTY. 


\ 


1. — Come & noto, le equazioni fondamentali della cascata nell’approssi- 
mazione B di Rossi e GREISEN sono: 


ön(E, t) ön(B,t) _ 
(Lo) ô 2. DE m 
co E E 
E'— E r / { r 
= fa, 1)R(E', EE) - an dE'— x(E, H fre, E') a db’ + 
E 0 
dH’ 
+2 [vor t)R(E, E') E’ 
E 
ö B,t ! ! n] E ! 
(1b) me JB: [zw ‚t)R(E', E) a dE' — Dy(E, t), 
E 
con 
4 { E' NOT: 
(10) Rw, B) =1—(«+ 5) 5 (1-5 ee 
dove 
(1d) & = 9 In (1837-12) ~ 1/40. 


Per x(H,t) (rispettivamente y(Z,t)) s’intende il numero medio di elettroni 
(rispettivamente fotoni) di energia compresa fra Æ e H+ dÆ, nell’intervallo 
di profondita che va da t a t+ dt, con ¢ misurato in unità di radiazione. 
B rappresenta la perdita di energia per ionizzazione per unità t. 
Le sezioni d’urto per produzione di coppie e di Bremsstrahlung sono quelle 
asintotiche usate da Rosst e GREISEN, benché espresse in forma più compatta. 
Applicando la trasformata di Laplace 


Titi) [exp [— Atf(R, t) dt 


0 
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alla (la) e alla (1b), si ottiene il seguente sistema: 


SLA(E, A) 


(2a) AL„(E, A) — B SE = r(BE,0)+ 
> r 2 r 
+ [enr B nam LAB, nf mn, EB’) a dB’ + 
E 0 ‘ 
+2 | LE", ARB, BY), 
E 
= E 
(20) (+ DEAR, 2) = | Dal BY, ARB, B) ay, AB" 


B 


avendo posto y(H,0) = 0 e chiamando x(E, 0) la sorgente elettronica. 
Per risolvere il sistema (2), si noti che, tenendo conto della (2b), l’ultimo 
termine della (2a) si puö scrivere: 


CO CO 
ear oN 2 ER NN 
(3a) 2 | LE’, DRE, EB) ar = Fp | L,(E', AR(E', B) Zr à 
E B 
dove 
D, ni / dn 
(3b) R(E', B) = | RB’, n)R(B, n) ays 
E 
e percid la (2a) diventa: 
à dL;(E À) D ae nl 7 Bis 2 
(4) AL„(E,A)—B 2. = n(E,0)+ [28 an, 2 Dr dE 
E 
E co 
I 7 ( E ! 2 ny D E'. B dh’ 
STB, 2) | rin, By „an + an L;(B", ARE, B) a : 
0 E 


La (4) @ un’equazione integro-differenziale nella sola funzione incognita 
L,(E, À), essendo note x(E,0), K(E', E) e R(H', E). 


2. — Una delle condizioni di limite cui deve soddisfare la Z,(E, À) è che 
sia (?): 


(5) lim L,(E, À) 0 almeno come 1/H**", 
ù E— co 


dove s & un numero complesso la cui parte reale & > 0. 
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Ritenendo questo limite, e sostituendo nella (4) a E la variabile x, inte- 
griamo ambo i membri di questa equazione su x da H a oo, e abbiamo, cam- 
biando opportunamente gli ordini di integrazione: 


(or era) — 0 0) dæ — À | date À) da — 


7’ 


co 
E' E' / À / 2 ' 2 E'— 2 


2 | ; da 
Le „(Ei jaw [mo Ci 


Nell’ultimo e nel penultimo integrale, dove sono stati scambiati gli ordini di 
integrazione, E’ varia da Æ à co mentre x varia da H a LH’. 

Poichè é lecito scambiare nei primi tre integrali la x con la ZH’, si ottiene 
infine: 


co co 


(Ta) BEE AN il x(E, 0) aE’ — [L.(B, AK(E',E)aE' |, 


(7) KE,D=i+ | run. x) = da — 


E'— x 2 À 
— Ja E'— x) DE dx D ET x) 
7 E 


La (7a) rappresenta un’equazione integrale di seconda specie del tipo di Volterra 


di cui la (7b) è il nueleo. 
La trasformata di Mellin applicata alla (7a) porge: 


(8) | ei LAB, 2) a 5) nee 
= NG) eo (elek) free ann ni - 


M,,(s +1, 0) | AT a B 
en N,(s + 1, À) %) K(Z", als): 
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Calcoliamo Vintegrale a terzo membro ponendo v = E/E’ e il valore di 
K(E', E) dato dalla (7b), cioè (v. Appendice I): 


1 1 


(9) forear = I IN ic 2) — 1] R(z)z dz — 
s EM 


0 
2 : ZN sit 2 : 
u ya  — — = = s+2 q 
A+ a | (5) R(z, w) a PIE R(w) aw}. 
0 


0 


Negli integrali a secondo membro si riconoscono rispettivamente i coefficienti 
A(s +1), B(s +1) e C(s +1) di Rosst e GREISEN. Percid, chiamando /,(s) 
e /,(s) le radici dell’equazione: 


(10) 2? + A[A(s) + D] + A(s)D — B(s)O(s) =, 


si ha: 


1 


dar for _ 1 Baht DA As +10) et 


srl AD sr joe 


0 
e la trasformata di Mellin della (7a) diventa: 


M,(s +1, 0 G(s +1,4 
(12) WAG: 2) = = a ) NAG Se il; À) it ’ 


che, tenendo conto delle diverse notazioni e della integrazione eseguita sulla (4), 
si riduce facilmente all’eq. (5) di Scott. 


3. - Un rapido esame del nucleo (7b) ci mostra che tale nucleo non & né 
simmetrico ne antisimmetrico. Possiamo dunque fin d’ora affermare che la (7a) 
non ammette autovalori e quindi non sara possibile ottenerne la soluzione con 
serie di autofunzioni, usando i noti metodi variazionali (°). 

EB noto perd che l’eq. (7a) ha una e una sola soluzione agevolmente ricava- 
bile con il metodo delle approssimazioni successive. Con tale metodo, ponendo 
al posto di L,(E', À) il primo termine a secondo membro come approssima- 
zione zero e, via via, i valori di L,(E, 4) che si ottengono successivamente a 
sinistra, si raggiunge rapidamente la soluzione: 


co co EK’ 
(13a) BLa(B, 2) = | (Bi, 08 + 3 (— a) | aH, oa [kun mar, 
n=1 n Ë 


E 


(6) G. Sansone: Eq. differenxiali nel campo reale. I (Bologna, 1948), p. 293 e segg. 
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con 
(136) K,(y, B) = | Katy, )K(@, B)de (n=2,8,4,.; Ky = K), 
E 


relazione, quest’ultima, che da l’n-simo nucleo iterato dal precedente K,_,(y, E). 
Nel caso di una sorgente elettronica monocromatica di energia Hy, cioè con 


(14) ra(E',0)= d(B, — E"), 
la (13a) si serivera: 
Ey 
co 1 n 
(15) BLa(Bos E; À) =1+ > (— 5) ET E) dy , 
n=1 
E 


e con una sorgente elettronica qualsiasi, per il principio di sovrapposizione, 
si avra dalla (15): 
> 
(16) L.(B, 4) = | (E', 0)Z.(E', B, 2) AB" 
Se vogliamo ora applicare 3a trasformata di Mellin alla (13a), tramite la 
(13b) e la (11), otteniamo (v. Appendice IT): 


: co ; s! : 
(17) N .(s, À) ak u (stn+1)! 


+ Ms +n-+ 1,0) G(s +1, A) G(s + 2, A)... G(s + n, À), 


che & l’eq. (10) di Scott, salvo le solite sostituzioni. 

La (15) non si presta all’antitrasformata di Laplace in quanto ogni suo 
n-simo termine ha 2n zeri distinti e un polo di ordine n nel punto À = — D 
e la medesima osservazione vale per l’eq. (10) di Scott. 

Questo Autore, perd, seguendo un metodo indicato da-SNYDER (!), per mezzo 
di prodotti infiniti delle funzioni G(s, A), con valori opportuni delle s, ha otte- 
nuto l’antitrasformata di Laplace, utilizzando gli infiniti poli, successivi ai 
2n zeri di ciascun termine, che a questo modo si sono venuti a creare. 


4. — Vediamo ora come la presente impostazione possa condurre alla solu- 
zione (14) di BHABHA e CHAKRABARTY. 
A tale scopo consideriamo una rappresentazione della ZL,(E, A) del tipo: 


co \n 1 k ey 
(18) L.(E, À) = > | a ial oe ds, 


n=0 
G 
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dove ( indica un percorso di integrazione che va da g—ico a 0 + ioo e cioè 

la retta parallela all’asse immaginario che passa per o>1, giusta la nostra 

condizione di limite (5) assunta all’inizio del n. 2. Evidentemente i coefficienti 

delle potenze di (—f/E) sono antitrasformate di Mellin da determinarsi. 
Tenendo conto della (18), la (7a) si pud anche serivere: 


(19) BLAB, 4) = 5 


2700 
C 


1 M„(s—1,0) ds 
s—1l [RE 


n=0 


co 

co B\" al REC: n) C/A \stn _ 

> | m) er as | a K(E', B) ab’; 
C E 


Ponendo come al solito v = E/E', l’ultimo integrale a destra, per la (11), 
si scrive: 


G(s + n—1, A) 
stn—1 


co 1 
1 
(20) E | za (*) = B [ext dv = E 
A ) J 


1! 


Dividendo ambo i membri per B, la (19) diventa allora: 


(STO) date 
Al Est | 


DTM 
ED LB, =a i 
C 


ee Bit 1 [ks,n) Gs+n—1, À) 
+3 ( 5) Qui | Es stn—1 GE 
6 


Sommando il primo termine a secondo membro con il primo termine. della 
serie, e scrivendo il resto di quest’ultima spostando n in n + 1, potremo espri- 
mere la (21) cosi: 


| ‘ree! | 
@2) DB ) =a = — {Mus —1, 0) — Ms, (6-1, A)+ 
6 


ds. 


> 


n=0 


Bi [= n +1) @(s+n, À) 
5) Zi, Es s+n 
é 


Tl secondo membro della (22) coincidera con il secondo membro della (18) se 


(23a) k(s, 0)G(s —1, À) = M;(s —1, 0) 
e 
(235) eh LE ee OR >) 


s+ 7” 
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Queste relazioni ricorrenti corrispondono, mutatis mutandis, alle relazioni (8) 
e (13) usate da BHABHA € ÜHAKRABARTY. 

Infatti, tenuto conto che la funzione w,(s,t) di questi Autori è un’oppor- 
tuna trasformata di Mellin della x(£, t), applicando alle loro eq. (8) e (13) la 
trasformata di Laplace, si ottiene rispettivamente, con i nostri simboli: 


(24a) N,,,(8, À) = 1/G(s —I, A) 
e 
(24D) Ny, (8, À) = Ny, _, (8 Hg + n—1, À) METER 


per il teorema del Faltung e per la (24a). 
Bastera ora porre in queste: 


k(s, n) (8 — 1)! 


= N =(, 12202 
M,(s—1, 0) (8+n—1)! (n ha 


(25) Ny, (8, À) == 
perchè le (24) si riducano alle nostre (23). 

Col solito procedimento ricorrente si trova dalle (23) il valore di k(s, n) 
che e; 


(s+n—l1)! M,(8 —1, 0) 


(26) Lan = (n = 0, 1, 2, «.). 
TI Gs + »—j, A) 


Ponendolo infine nel secondo membro della (18), otteniamo, per il caso di una 
sorgente monocromatica (M,(s —1, 0) = Hi-*): 


27 1 7 = pi 1 "fE,\(8 + n—1)! ds 
re UNS > (— Fa mix] (a) reger . 
À à [6] | [ G(s-+--n —j, A) 
f=1 


La (27) è la trasformata di Laplace della soluzione (14) data da BHABHA e 
ÜHAKRABARTY. 

Resta cosi dimostrato che la soluzione di questi Autori è una rappresenta- 
zione in serie di potenze di (— f/f) della nostra soluzione (15), di cui l’eq. (10) 
di Scott @ la relativa trasformata di Mellin. 


E mio piacere ringraziare qui il prof. Puppı per il suo cortese interessa- 
mento e il prof. CLEMENTEL per le utili discussioni nello svolgimento del pre- 
sente lavoro. 


J 
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APPENDICE I 


a) I coefficienti di Rossi e Greisen della (9) si trovano nel modo seguente 


@ = aL’): 


1 1 uf 
[rat | mr ET — R)(1 — 2) de | = 
; 2 


0 


8 
z)s+1 — 1]R(z) — AG sasha 
aos = se 
avendo cambiato l’ordine di integrazione e sostituito 1—z con z nel secondo 


integrale. 
b) Ponendo al posto di R(e ) il suo valore dato dalla (3b) e cambiando 


due volte Vordine di integrazione, si ha (w = n/E'): 


1 3 
2 [var | Re) dz = 
0 v 
1 


1 \ 
= sone [os a)" R(z, w) af) = = Ah ae DE Be ae À 
8 +1 Ww w Fe 


0 0 
poiche 
1 1 
2 \st1 2) 
B(8 + 1) = —) Æ(e, w)d A e O(s + 1) = | w**R(w) du. 
w w 
0 0 


APPENDICE II 


Dividendo ambo i membri della (13a) per 6 e applicandovi la tr ‘asformata 
di Mellin, cambiando, nell’n-simo termine della serie, l’ordine di integrazione 


due volte, si trova: 


[6 ) 2.08, na) Sone = 


co Bi y 
_ Me +1,0) ,1S/. 1" x a (=) a By oe 
aaron wall a rw oar fa (ro » : 
0 0 0 


— 
— 
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Moltiplichiamo ora ambo i membri della (135) per (Z/P)’d(E/ß) e inte- 
griamo da 0 a y; cambiando a secondo membro l’ordine di integrazione una 
volta soltanto: 


uy Y 


Es a @(: inal srl a 
(5) TR To Hy Er AE 6) Kt 


B B s +1 B Pe 


avendo moltiplicato e diviso K(z, H) per (4/B)°*1 e usando la (11). 
Questa relazione, per ricorrenza, ci da subito: 


Y 


ENT Ve s! DONE ae 
fo) K, (y, E) B = ji (3) Ga il G(s + 1, A)G@(s + 2,4)... G(s + n, À), 


0 


che, posta nella (-), la riduce alla (17) del testo. 


SUMMARY 


By using Laplace transformations, the diffusion equations of the shower theory are 
reduced to a Volterra integral equation, a solution of which is obtained by using the 
method of successive approximations. It is shown that this solution corresponds to 
Scott’s one. Further it is shown that the same equation leads to the solution of 
BHABHA and CHAKRABARTY. 
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Sulla teoria evaporativa delle disintegrazioni nucleari. 


F. FERRARI 
Istituto di Fisica dell’ Universita, Centro di studio degli Toni veloci del C.N.R. - Padova 


C. VILLI 


Istituto di Fisica dell’ Universita - Trieste 


(ricevuto il 7 Aprile 1952) 


Riassunto. — Il modello a gas di Fermi é stato esteso alle alte tempera- 
ture nucleari mediante un procedimento matematico non soggetto alla 
nota condizione restrittiva t/y < 1. I processi evaporativi per varie energie 
di eccitazione del nucleo sono stati studiati tenendo conto dell’effetto 
delle interazioni nucleari, della diminuzione della temperatura nucleare 
e della variazione, durante il decorso del processo, del numero di massa. 
Si discutono le previsioni teoriche che risultano in soddisfacente accordo 
con i dati sperimentali. 


1. — Introduzione. 


Lo studio delle disintegrazioni nucleari viene attualmente compiuto in gran 
parte mediante l’esame di stelle nelle emulsioni fotografiche. I numerosi risul- 
tati sperimentali giä raccolti con questa teenica hanno dato lo spunto a diverse 
considerazioni teoriche intese ad interpretarne le modalita più importanti. 

Limitandoci qui al solo caso di disintegrazioni provocate in nuclei pesanti 
dell’emulsione (Ag, Br), il quadro teorico puö essere riassunto nel modo se- 
guente. Il nucleone primario incidente, ionizzante 0 no, interagisce in un primo 
tempo con un solo nucleone del nucleo e puö dare luogo sia ad un urto ela- 
stico sia ‚se la sua energia & sufficientemente elevata, ad uno anelastico. In 
questo caso vengono emessi un certo numero di mesoni, che appaiono in genere 
come rami di sciame, e nucleoni di rinculo e di knock-on espulsi dal nucleo 
per urti secondari dai mesoni e dalla coppia primaria di nucleoni. Nel caso, 


488 F. FERRARI € C. VILLI 


invece, dell’urto anelastico si hanno ancora generalmente nucleoni di knock-on, 
ma mancano i rami di sciame, 

Con cid ha termine la prima fase, quasi immediata del processo; succes- 
sivamente> ha inizio, un processo molto più lento dovuto al fatto che parte 
dell’energia del primario viene ceduta per urti larghi ed & rapidamente distri- 
buita a tutti i nucleoni del nucleo. Il conseguente stato di eceitazione nucleare 
diffusa provoca l’emissione di rami di bassa energia (tracce nere) secondo il 
modello evaporativo proposto da WEISSKOPF (1). 

Recenti tentativi (FERMI (?)) sembrano pure aprire una via promettente 
per Vinterpretazione delle collisioni anelastiche e rendere conto dell’anda- 
mento delle particelle di sciame. Per quanto riguarda invece la zona delle 
energie intermedie (rami grigi) si hanno per ora solo indicazioni qualitative 
e vi sono anzi a tale proposito divergenze di opinione che fanno apparire la 
loro interpretazione ancora piuttosto problematica. 

Un primo modello, proposto da GOLDBERGER (*), considera i rami grigi 
come essenzialmente dovuti a protoni di rinculo e di knock-on nell’urto diretto 
del primario e negli urti immediatamente successivi a questo. Secondo questo 
schema, Venergia del primario E, viene spesa in massima parte oltre che nella 
creazione dello sciame di mesoni, nell’emissione di rami grigi di knock-on in 
modo che Venergia di eccitazione del nucleo Hy, responsabile del fenomeno 
di evaporazione, à solo una piccola frazione di Æ,. (Anche quando non vi 
sono particelle di sciame essa non supererebbe i 100 MeV). 

Il modello opposto, prospettato da CAMERINT e col). (*5) sulla base dei 
loro: dati sperimentali relativi a stelle di cospicua energia, pure non esclu- 
dendo che parte dei rami grigi possa effettivamente corrispondere a nucleoni 
di rinculo e di knock-on, considera la maggior parte di essi come dovuta ad 
un prolungamento dello spettro di evaporazione verso energie più alte di quanto 
la solita teoria lasci prevedere. Segue allora che, in assenza di rami di sciame, 
Venergia di eccitazione #, potrebbe essere prossima all’energia ceduta al nucleo 
dal primario EP, e, secondo il gruppo di Bristol, raggiungere valori intorno 
a 800-1000 MeV. 

Senza voler affatto discutere quale dei due suddetti modelli sia da consi- 


) V. Werisskorr: Phys. Rev., 52, 295 (1937). 
) E. Fermi: Progr. Theor. Phys., 5, 570 (1950). 
3) M. L. GOLDBERGER: Phys. Rev., 74, 1269 (1948). 
) U. Camprint, T. Coor, P. H. Fowrer, W. O. Lock, H. MurRHEAD, e N. ToBın: 
Phil. Mag., 40, 1073 (1949); U. Camprint, P. H. Fowver, V. O. Lock e H. MuIRHRAD: 
Phil. Mag., 41, 413 (1949); U. CAMERINI, J. H. Davies, P. H. FOwLEr, C. FRANZINETTI, 
H. MuviıkneAnD, W. O. Lock, O. H. Perkins 6 G. YEKUTIELY: Phil. Mag., 42, 1241 
(1951). 

(5) U. CAMmERINI, V. O. Lock e D. H. Prrkıns: Progr. in Cosmic Ray Phys. (Pub.) 
(1952). 
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siderarsi più aderente alla realtà fisica, abbiamo carcato di modificare l’usuale 
teoria della evaporazione per stabilire se essa possa rendere conto dell’emis- 
sione di almeno parte dei rami grigi, in modo da interpretare lo schema pro- 
posto dal gruppo di Bristol meglio di quanto sia stato sinora possibile. 

Com’& noto, l’andamento dello spettro evaporativo delle particelle emesse 
dal nucleo dipende essenzialmente dalla relazione postulata tra energia di ecci- 
tazione, temperatura nucleare ed entropia. Tutte le tecrie sinora svolte, salvo 
poche eccezioni basate sul modello a goccia, si sono fondate a tale scopo sulle 
formule della statistica di Fermi nel caso limite di un gas totalmente dege- 
nere, ciod per 7 cosi prossimo a zero che la relazione tra E, e T possa scriversi 
7, = b-T* con b costante. 

Anche i miglioramenti introdotti da vari Autori (°?) alla primitiva teoria 
© consistenti nel tener conto della diminuzione della temperatura nucleare nel 
corso delle successive emissioni di particelle, nel considerare le forze di inte- 
razione fra i nucleoni, ecc., poggiavano sempre su questa medesima ipotesi. 
Tuttavia, ei si pud facilmente convincere che, nel caso di disintegrazioni con 
energia di eccitazione rilevante, la corrispondente tempe satura nucleare & di 
un ordine di grandezza tale da rendere molto inesatta la relazione Ey = b-T? 
e quindi non appare lecito per queste energie applicare le formule che si rife- 
riscono ad un gas totalmente degenere, ma sembra consigliabile sostituirle 
con altre valide per qualunque temperatura. 

Con questa rettifica siamo riusciti a modificare verso le alte energie la 
forma dello spettro di evaporazione con un risultato che sembra, rispetto a 
quelli ottenuti da altre teorie, in migliore accordo con i recenti dati sperimen- 
tali, semprechè si accetti di interpretare parte delle tracce grige come un 
annesso del processo evaporativo. 


2. — Modello a gas di Fermi. 


Com’ noto, l’equazione di stato di un gas totalmente degenere obbediente 
alla statistica di Fermi si ottiene eliminando fra le due relazioni (°) 


Co 
(la) N= 2 c| E®f(H)dR, 
oat: 
3 N : 13/2 nl ni 
(1b) Utot = 5 ( | Esf(E)dE, 
ö 


(6) E. Bacar: Cosmic Radiation (Heisemberg, 1946), pag. 128. 
(?) K. J. LæCoureur: Proc. Roy. Soe., 63, 259 (1950). 
(8) H. A. Berne: Rev. Mod. Phys., 9, 83 (1937). 
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il parametro A= exp [—p/r] che compare nella funzione di distribuzione f(E) = 
= [Aexp [H/t] + 1]>. A tale scopo, si procede generalmente come segue: 
integrando per parti la relazione 


(2) 1, = | om ar, G12) 


in cui possono compendiarsi le (1; a,b) [®,(#) = E!:; ®,(E) = #2], si ot- 
tiene, dopo aver posto ®,(#) = — dF;/dE ed eseguito il cambiamento di varia- 
bile «= (E— y)/t 


(3) 


Nell’ipotesi che il gas sia fortemente degenere (A= exp [—yp/t] <1) la funzione 
Fi(y + tx) che compare nella (3) pud essere sviluppata in serie in un intorno 
di æ = 0 e, ponendo senza sensibile errore w/t — oo si ricava che l’energia di 
eccitazione Z, è legata alla temperatura nucleare 7 dalla relazione H, = b-7?, 
in cui b è una costante (°). 

Per ottenere una relazione tra EZ, e t che sia valida per ogni temperatura 
nucleare, & necessario ricavare l’equazione di stato del gas nucleonico indi- 
pendentemente dalla condizione restrittiva t/py<1. A tal fine, posto y = t/y, 
le (1; a, b), tenuto conto della (3), possono scriversi nella forma: 


al e)3/2ex 
(da) N do = Opal) 
—ı]u 
ne } 
6 (A + ya)5l2ex Gas 
(4b) Utot = 5 Cys | à ae da = - Cysy(y) . 
—ıly 


Fissato il valore di y, la (4a) consente di calcolare il valore del parametro w 
in funzione di y (*) (2 = N = A/2) 


20 


A \2i8 
(4e) ply) = | Bey), 


(2) J. E. Mayer e M. G. MAYER: Statistical Mechanics (Oxford, 1938), pag. 374. 
(*) La p(y) definita dalla (4c) à una funzione decrescente con la temperatura +: 
nella zona di Fermi (y < 0, 1) tale diminuzione avviene secondo la legge 


Y = poll — (22/12) (t/po)?] : 
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che, sostituita nella (4b), da: 


A \518 
(4d) Une (5) we 8(y)v(Yy) - 


2 


Dalla (45), tenendo conto che u(0) = »(0) =1, si ricava per t= 0 (y= 0) 
Venergia minima del gas di Fermi 


60e 
(de) Tg WR 
essendo per la (4c) 
A \2/3 
4(f) yO) == (55) 


che esprime l’energia di Fermi caratteristica degli insiemi che obbediscono 
al principio di esclusione di Pauli. 
L’energia di eccitazione 1isulta quindi: 
= 60)/( A\*s 
(5) Bo = Ue) 00 =" (56) ery) — 4. 


5 \2C 


La (4e) definisce la massima energia cinetica sia dei neutroni che dei pro- 
toni, misurata a partire dal fondo della buca di potenziale. In virtü della 
legge R — 7,A¥%, in cui R é il raggio del nucleo di numero di massa Ase 
r, = 1,5-10-13 cm è il raggio di un nucleone, si ha che la costante C varia 
linearmente con il numero di massa e y, assume, tanto per i protoni quanto 
per i neutroni, il valore 21 MeV. Tenendo invece conto dell’eccesso neutro- 
nico si trova che l’energia di Fermi per i neutroni e, rispettivamente, per i 
protoni è: 


OTA BR DA 2/3 2 7 2/3 
(4h) Yon — =) Wo 3 Vor — | 7 | Vo - 


3. — Interazione dei nucleoni. 


La relazione (5) & valida nell’ipotesi di un gas ideale. Per descrivere invece 
in modo più adeguato il comportamento del gas nucleonico terremo conto 
del fatto che le Z particelle cariche e gli À — Z neutroni del nucleo di numero 
di massa À interagiscono fra di loro. Di queste interazioni trascureremo quelle 
di natura coulombiana dovute alla repulsione elettrostatica fra i protoni, in 
quanto la correzione sull’energia totale che esse comportano non é molto sen- 
sibile neppure per i nuclei pesanti (12). 


(10) R. D. Present: Phys. Rev., 60, 28 (1941); E. FEENBERG: Phys. Rev., 60, 
593 (1941). 
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L’interazione nucleare propriamente detta potrebbe essere introdotta come 
una perturbazione e trattata con gli usuali metodi variazionali in cui il raggio 
nucleare viene considerato come una quantita variabile. Noi tuttavia ci rife- 
riremo ad un eriterio diverso, sebbene sostanzialmente analogo e più adeguato 
al modello a gas di Fermi. 

Se si tiene conto delle interazioni nucleari, lo stato di estrema degenera- 
zione del gas non pud essere descritto dicendo che tutti i livelli più bassi sono 
occupati, in quanto, a causa delle transizioni che esse producono fra i livelli 


10 


AS a. 
at 
~ 


——— Distribuzione di Fermi per T= 0 

— - — - Distribuzione perturbata dalle in- 
terazioni nucleari. 

Karmann Distribuzione di Fermi corrispon- 
dente alla «temperatura equiva- 
lente ». 


Fig. 1. — Distribuzione dei momenti in un nucleo pesante standard. 


oceupati e quelli non occupati, la distribuzione di Fermi viene distorta in un 
intorno di y, nel senso che nello spazio dei momenti la densitä tende a dimi- 
nuire per p< yi/2m e presenta una coda per p> yi/2m (fig. 1) (1). Si pud 
quindi pensare che l’effetto delle interazioni nucleari sia dovuto ad una specie 
di «agitazione termica » della materia nucleare; si prospetta in tale modo la 
possibilità di individuare per ogni temperatura effettiva 7 del nucleo alla quale, 
a causa delle interazioni nucleari, la distribuzione di Fermi & perturbata, una 
temperatura equivalente r che caratterizzerebbe una distribuzione di Fermi 
relativa a particelle che non avessero interazione. 

Eseguendo tutti i ealcoli tenendo conto esclusivamente della temperatura 
equivalente e risalendo poi a quella effettiva secondo una determinata legge, 
si ha la possibilità di ineludere le interazioni nucleari nello studio del com- 
portamento del gas senza introdurre aleun criterio fisico 0 accorgimento di 
‘alcolo estraneo alla statistica di Fermi. 

Per far cid & necessario stabilire anzitutto il legame che intercede fra la 
temperatura effettiva del nucleo e quella equivalente. A tale scopo, riguarde- 
remo per le ragioni dette le relazioni (4; a, b) come dipendenti dalla tempera- 
tura equivalente + ed essendo la y una funzione decrescente della temperatura 
equivalente + alla quale la distribuzione di Fermi & realizzata, faremo l’ipo- 


(41) L. RoseNreLp: Nuclear Forces (North-Holland Pub. Co., 1948). 
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tesi che tra le due temperature 7 e + sussista la relazione 
(6) t= T + ay(t), 


dove, secondo WATANABE (#338), a = 0,29. 
Allora, la (4e), interpretata come funzione di +, diventa per T= 0 


L A \2/8 

(Ta) WT) = oh uS(a), 
cioè 

(Tb) To = Ay(To) ; 


e quindi la relazione che sostituisce la (4e) risulta: 


of h 6C/ A \5S . 
(Te) Il), Bea) 5 
o \20)] — 
la differenza fra la (Te) e la (de) d& Venergia latente di eccitazione nucleare 
dovuta alla equivalenza termica delle interazioni fra i nucleoni: 


GC / A \sis 
U'— ar ( | {uSS(a)r(a) — 1 Le 


2C 


EB di immediata verifica che per ogni valore del parametro y, cui cor- 
risponde una determinata temperatura equivalente, si trova la temperatura 
effettiva del nucleo tramite la relazione: 


A 


2/5 
—2/3/ 
20 | u (y) : 


(7d) T = (y— à) | 


in cui y >a. Inoltre, fra temperatura equivalente, quella effettiva e il grado 
di degenerazione del gas esiste la relazione: 


(Te) t= Til + a log Ay. 


Tenendo conto delle interazioni nucleari, la relazione (5) diviene per i neu- 
troni e, rispettivamente, per le particelle eariche: 


3C (2A — 2Z\58- $ 

(Sa) ze = ( en) Lee S(y)o(y) — u SS(a)»(a)| , 
30 /2Z\8S 5 IST 

(Sb) Lo > T (=) (ud) — ee SS(a)r(a)| , 


(4?) S. WATANABE: Zeits. f. Phys. 113, 482 (1939). 
(3) H. Eurer: Zeits. f. Physik, 105, 553 (1937). 
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e quindi Venergia totale di eccitazione del gas nucleonico risulta 
(10) By = Kart) = List) + H(t). 


La (9) @ stata calcolata numericamente con grande approssimazione per 
diversi valori del numero di massa A profittando del fatto che i contributi 
agli integrali ply) e v(y) sono 
del tutto trascurabili per # > 30 
ed y < 10%, a causa del massimo 
acuto della funzione e”(l + ¢)-* 


0 


> 
a 
zZ 
€ 
€ 
uy 
‘o 
2 


in corrispondenza di #0. 

Il risultato dei calcoli € rap- 
presentato in fig. 2 in cui l’ener- 
gia di eccitazione HM, & stata 
posta in funzione della tempera- 
tura effettiva del nucleo tramite 
la (7d). Com’era da attendersi, a 
causa del fatto che le interazioni 


Fig, 2. - Curve energia-temperatura per nuclei 
di differente numero di massa. Le curve trat- 
teggiate rappresentano la variazione di energia  nucleari tendono a ridurre lo 
in funzione della temperatura nucleare durante stato di degenerazione del gas 
il decorso del processo evaporativo iniziatosi con nucleonico, le curve energia tem- 

E, = 800 MeV e E, 200 MeV. 


peratura tendono a raccordarsi 


nella zona delle alte temperature 
con quelle calcolabili per gas di media degenerazione con il metodo di Dar- 
win-Fowler, 
In base alla fig. 2 Ventropia del gas nucleonico à stata calcolata tramite 
la relazione 


E 


(10) S(K) | na | 


0 


Il risultato dei calcoli & rappresentato in fig. 3. 


4. — Evaporazioni nucleari. 


Secondo la teoria dell’evaporazione di Weisskopf (), la probabilità che il 
nucleo eccitato al livello Hy, emetta una particella di massa m con energia cine- 
tica compresa tra ee + de, à: 


Smym 
(11) Wie) de = o, mh exp [SCH, U — e) — S(H,)lede, 
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in cui y rappresenta il fattore statistico di spin e U Venergia media di legame 
della particella del nucleo, Nella (11) abbiamo supposto costante la sezione 
d’urto o, del processo inverso di quello deseritto dalla (11) e eid non comporta 
un serio errore sin tanto che Venergia della particella evaporata si aggira sui 
25-30 MeV (rami neri). 
L’energia einetiea media risulta 
liy-U M—U 


(12) e | eW(e) de /| W(e)de. 


0 


Per basse energie la maggior parte delle particelle emesse da un nucleo 
dopo la collisione sono relativamente lente; la distribuzione energetica delle 
particelle cariche ha un massimo che corrisponde approssimativamente allal- 
tezza della barriera di potenziale, mentre le particelle neutre hanno una distri 
buzione maxwelliana con il magsimo alla temperatura del nucleo, 

La (11) à stata ealeolata tramite i valori ricavati dalla (9) @ dalla (10). 
A causa del massimo acuto presentato in uno stretto intorno di 4, - U dalle 
funzioni integrande e exp [S(M,— U--e)| e et exp [SUM U K)] che com- 
paiono nella (12), à facile verificare che Venergia einetien media della s-sima 
particella evaporata à inversamente proporzionale al eoeflieiente angolare 
k( Hy — U) della curva energia-entropia caleolato in corrispondenza dell’energim 
posseduta dal nucleo che evapora tale particella: 


(13a) &, = 2/k(,.,— VU), 
6, quindi, tenendo conto che dS/d/7 1/7’ = k, si ha 
(135) @, = 2T(D, nu — VU). 
Dette 
iy by 


Z ‘ N 
: | Wie) de, Py N | W y(e) de, 


iy, i 


(14) P, 


la probabilità totale di emissione di una particella carica 6, rispettivar 
mente, di una particella neutra, Venergia media di legame che compare nolla 
(12) & data da; 


(15) U U(k,) U, | = Uy PM) U, | y( My) Uys 


(4) N. Pac: Proc, Phys. Soe., 68, 250 (1950). 
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in cui U,e U, sono state ricavate dai dati di WEIZSÄCKER 


U, = 3) MeVe: = Ui 1452 Mee BU 219) 3 MeN. 

e tenendo conto del fatto sperimentale che il rapporto fra i protoni e le altre 
particelle di carica unitaria (deutoni, tritoni) e 55/79 = 0,696 (14). Abbiamo 
trascurato l’emissione di particelle alfa in quanto, se non si considera l’effetto 
tunnel, il rapporto P,/Py è molto piccolo rispetto a quello relativo a qual- 
siasi altra particella carica; questa semplificazione, che altera il rapporto tra 
i singoli tipi di particelle cariche, non incide che in maniera poco rilevante 
sul numero totale di particelle ionizzanti. 

A causa del comportamento delle funzioni integrande che compaiono nella 
(14) e per ragioni analoghe a quelle precedentemente esposte, le P; e Py 
risultano inversamente proporzionali al quadrato del coefficiente angolare della 
curva energia-entropia calcolato nei punti #,— U, e H,— Uy. Tenendo conto 
della relazione 1/7 = k, l’energia media di legame della s-sima particella eva- 
porata si ricava allora in modo ricorrente dalla (15) e si trova: 


= A(s) TE, — U,)U, + pls) T?(B,_,— Uy) Uy 
A(s)T?(B.1 — Us) + ws) TE, — UN) 


(16) 


U (Eo) > 


’ 


in cui 


| Als) = Z—(s—1)D(E, ;); 
| w(s) = N= (8 —1)x(E,_,). 


Il calcolo numerico è stato 
eseguito seguendo le singole fasi 
del processo evaporativo. Fis- 
sato il numero di massa A e 


© 
>> 


(17) 


10° E, in MeV 


ie) 1 2 3 4 5 6 27. 8 9 


Fig. 3. — Curve entropia-energia per nuclei di 

differente numero di massa. Le curve tratteg- 

giate rappresentano la variazione di entropia 

in funzione della energia di eccitazione del nu- 

cleo residuo durante il decorso del processo 

evaporativo iniziatosi per H, = 800 MeV e 
E, = 200 MeV. 


l’energia di eccitazione E, del 
nucleo, l’energia cinetica delle 
successive particelle emesse à 
stata valutata tramite la (12) e 
la (16). Nelle fig. 2 e 3 le curve 
tratteggiate rappresentano le va- 
riazioni dell’energia totale di ec- 
citazione del nucleo e dell’en- 
tropia durante il processo evapo- 
rativo iniziato con energia di ec- 


citazione di Æ,—800 MeV e E,= 200 MeV. I calcoli sono stati eseguiti anche 
per valori di H,=:400 MeV e E,=600 MeV. Nella fig. 4 sono rappresentate 
le energie cinetiche delle particelle evaporate in funzione dell’energia di ec- 
citazione. 
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Le successioni numeriche delle perdite di energia del nucleo per particella 
evaporata 7, = &+ U, sono delle successioni aritmetiche di ordine superiore. 
Tenendo conto delle approssimazioni introdotte dal calcolo, esse possono rite- 


nersi individuate da progressioni 
aritmetiche del terzo ordine; cid 
sarebbe ancora più esatto qualora 
si fosse tenuto conto che il rap- 
porto tra energia cinetica media 
acquistata dalla particella evapo- 
rata e l’energia totale a disposizione 
dipende oltre che dal modello as- 
sunto anche dalla forma o(E,, €) 
della sezione d’urto del processo 
inverso di quello descritto dalla (11), 
che noi abbiamo supposto costante, 
mentre invece & una funzione de- 
crescente, al crescere dell’energia &, 
come et, esclusa la zona delle basse 


30 


20 


| 


L 


energia cinetica in MeV 


10? E in MeV 


energie dove interviene il principio 
di esclusione. 
Nella tabella I sono indicate per 


Fig. 4. — Energia einetica media delle parti- 
celle espulse in funzione dell’energia di ecci- 
tazione per processi evaporativi iniziatisi con 
E, = 800, 600, 400, 200 MeV. Le linee trat- 


le diverse energie di eccitazione le 
perdite energetiche del nucleo per 
evaporazione della prima particella 
e le differenze prime, secondo, terze e 


teggiate rappresentano le isoterme relative 


alla temperatura vera del nucleo. 


quarte delle relative successioni nu- 


meriche. Le differenze quarte sono state considerate nulle in quanto il loro 
ordine di grandezza é inapprezzabile nei nostri calcoli. 


TABELLA I. 


mee N Anı At | An, | Atm, | A B | C | D 

| | | | | | 
800 | 40,4980 —0,4177 —0,9158 |—0,0049 | ~ 0 |—0,0049 | — 0,0338 |—4,8767 | 976,8674 
600 | 36.3570 — 0,4798  — 0,0305 | — 0,0049 0 —0,0049 0,0926 | — 5,4555 878,1110 
400 | 31,5250 — 0,6312 | — 0,0550 |— 0,0049 ~ 0 — 0,0049 —0,1906 6,9731 | 763,7686 | 
200 | 25.5001 —1,0053 | —0,0844 | — 0,0049 | — 0 —0,0049 |— 0,3082 | 11,1047 | 632,4278 


Poichè la somma S delle perdite energetiche che subisce il nucleo nel pro- 


cesso evaporativo completo è: 


N 
S=) 7. = Bo Um 
a 
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si ottiene immediatamente una formula di somma per una progressione aritme- 
tica del terzo ordine 


N N N N 
(18) Ey— U lm + (9) 4m + (5) aim + (7) Am. 
e quindi 
(19) FN) AN4 + BN3 + CN? + DN — 24(E, — U) = 0. 


I coefficienti A, B, OC, D si calcolano mediante le seguenti relazioni 


m 
Ann, = Nm ET MN m 4p (een Tee = (= 1)", ’ 
(204) 


m 
| Nm un + m An, IF (Jar 5 cise ar Am, ’ 
e si trova 
[4 = Ay, , CO = 12An,—12A%, + 11A%, , 
(205) 
| B — 4A%,— 6A’, , D = 24n, —12An, + 8A’, — 6A%, . 

La soluzione dell’equazione di quarto grado (19) @ alquanto laboriosa e 
dalla formula risolutiva non à possibile ottenere una relazione semplice e suffi- 
cientemente approssimata tra il numero totale di particelle evaporate e l’energia 
di eccitazione del nucleo. Per trovare una relazione non molto complicata 
abbiamo risolto la (19) mediante il metodo di Newton-Fourier, assumendo come 
radice approssimata per eccesso n — VE, (*). Siecome le condizioni per la 
applicabilità di detto metodo sono verificate, una radice della (19) è data da: 


La formula (21) ha perd l’inconveniente di richiedere la conoscenza delle 
A, B, ©, D che dipendono dall’energia di eccitazione termica del nucleo, come 
si pud rilevare dalla tabella I. Per questa ragione abbiamo cercato di modifi- 
care la (21) confrontandola con la soluzione rigorosa della (19); abbiamo con- 
statato che il rapporto F(VEE,)/F'(VE,) pud esprimersi, a meno di qualche 
unità sulla seconda cifra decimale, tramite la relazione: 


F(VE,) UE? —0,26 By 


22 E)o - ——— , 
(22) ÎCE) FP’ (VE,) Da 7 Hie 


(*) Una relazione di questo tipo pud ricavarsi dalla teoria di Weisskopf per nuclei 
pesanti aventi energia di eceitazione My, = 100 MeV (W. HeiItLer: Phys. Rev, 54, 
873 (1938)). 


SULLA TEORIA EVAPORATIVA DELLE DISINTEGRAZIONI NUCLEARI 499 


Tenendo conto della (14), il rapporto fra il numero delle particelle cariche 
e il numero totale di particelle evaporate in funzione dell’energia di eccita- 
zione E, & dato da 


A(s)T?(E.-,1 — Ui) 
A(s)T*(H,_, — U) + ws) TXE..— Uy) 


(23) DE.) = 
Segue che il numero delle particelle ionizzanti 
evaporate da un nucleo eccitato al livello E, è: 
(24) N (Ey) = D(B)VE, — P(£o)f (Eo) ; 


da cui si ottiene 


(25) qE) N; + PEN; — 2, = 9, 
essendo 
| Es 
Q(Eo) = => 5 ; 
a ® an ed (Bu)] 
Air eee 


(By) (Eo — PIE] 


Le funzioni f(E,), q(£o), € P(Eo) sono rappre- 
sentate in fig. 5. Fig. 5. 


5. — Considerazioni sulla legge E, = b*T?. 


Potrebbe essere utile, allo scopo di confrontare questa teoria con quella 
di Weisskopf (15), stabilire a quali conseguenze conduca l’espressione gene- 
rale (9) quando si preeisi la dipendenza tra energia di eccitazione Hy e tempe- 
ratura nucleare T mediante una legge del tipo: 


(27a) Hy —0T". 
Se si interpreta il nucleo come gas ideale totalmente degenere (7— 0), la (27a) 


con n = 2 & conseguenza diretta della statistica di Fermi e la costante b, legata 
alla struttura del nucleo dalla relazione 

4 8/8 V2/3m 3 Nea 
(27b) D = ; [((A—= Z)48 + ZU], 


Ta E 


assume, per À — 100 e Z — 44, il valore 11 MeV-!, Tuttavia, questo valore 
(5) V. F. Weısskorr: Lecture Series in Nuclear Physics (U.S. Govt. Printing Office, 
Washington, D.C., 1947). 
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non rende conto in modo adeguato dei dati sperimentali i quali impongono, 
secondo WEISSKOPF, per b un valore compreso tra 5 e 20 MeV-1, 

Nel nostro caso, tenendo conto che l’effetto delle interazioni nucleari con- 
duce ad una diminuzione dello stato di degenerazione del gas (1%) nella (27a) 
tanto b quanto n devono intendersi cariabili con H,. Per semplicità abbiamo 
considerato le due ipotesi estreme: 


MCE )E b= cost. = nh). 


In entrambi i casi si trova che b(Æ,) e n(E,) sono funzioni decrescenti del- 
Venergia di eccitazione. Da cid segue che nelle formule da noi ricavate la 
temperatura effettiva del nucleo appare essere minore nella zona di Fermi 
di quella che si avrebbe applicando le formule di Weisskopf, fatto questo 
che è in accordo con le osservazioni sperimentali di GUGELOT (17). Piü preci- 
samente, assumendo una legge del tipo (27a) si trova che per n = 2 e A = 100 
la funzione b(E,) varia da 9,9 MeV-1 (T= 6) a 3,6 MeV (T= 14,6 MeV). 

Ci sembra ancora interessante osservare che, tenendo conto della varia- 
bilita della massa durante il decorso del processo evaporativo iniziatosi, per 
esempio, con una energia di eccitazione termica 2, = 800 MeV, b(E,) varia 
tra 3,6 MeV-1 (T= 14,6 MeV) e 6,7 MeV-1 (T= 0). La dipendenza di b(#,) dal 
numero di Massa cosi calcolata conduce, anche in questo caso, a valori che 
nella zona di Fermi sono inclusi negli errori sperimentali di Gugelot. 

La seconda ipotesi b = cost. = 11 MeV !e n = n(E,) porta per un nucleo 
di massa 100 a valori di n variabili tra 1,96 (T— 0) e 1,66 (T= 14,6 MeV) (*). 


6. — Confronto con i dati sperimentali e conelusioni. 


La relazione (25), ricavata teoricamente, è stata confrontata con la for- 
mula empirica di Bristol (1) 


(28) Q'Ni+pNi—E =0, (qe 4 p 51) 


ottenuta in base all’analisi di una selezione di stelle prevalentemente di 
tipo 0,. Il numero medio di rami ionizzanti di natura evaporativa in funzione 


(46) J. BARDEEN: Phys. Rev., 51, 799 (1937). 

(7) P. C. GUGELOT: Phys. Rev, 81, 51 (1951). 

(*) Le funzioni a(H,) e n(E,) decrescenti per valori dell’energia di eccitazione cor- 
rispondenti a y >0,1 (n. 2) sono state approssimate con una costante nella zona di 
Fermi (y <0,1). 

(48) U. CAMERINI, J. H. Davies, P. H. FOwLEr, C. FRANzINETTI, W. O. Lock, 
O. H. PERKINS e G. YEKUTIELY: Phil. Mag., 42, 1241 (1951). 


RE» 
| 
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dell’energia di eccitazione termica del nucleo è in soddisfacente accordo con 
i dati sperimentali (fig. 6). 


a Dati sperimentali di PERKINS 16 


A » » di Bristol 
———— Curva teorica di LE COUTEUR 
= — » » di BAUMHOFF 
u » » di FUJIMOTO 
oe Znormular (29) 10 


numero di rami ionizzanti n; 


Fig. 6. — Curve teoriche (Le Cou- 
TEUR (7), BAUMHOFF (?4), FUJIMOTO (??)) 
e sperimentali (CAMERINI (17), PER- 
KINS (23)) del numero di ramiionizzanti 
in funzione dell’energia di eccitazione. 2 
La curva teorica di Baumhoff é solo 10° E, in MeV 
indicativa. 5 ae Soa. er OUR 


Nella fig. 7 (a, b, e) sono confrontati con i nostri (curva 2) gli spettri calcolati 
da LE COUTEUR (curva 1) delle particelle espulse dal nucleo avente energia 
di eccitazione H, = 250, 450, 650 MeV. Com’era da attendersi, le distribu- 
zioni calcolate in base ai criteri sui quali à fondata la presente trattazione si 
discostano sensibilmente da quelle previste da LE COUTEUR e mostrano come 
Vaumento dell’energia cinetica media dia luogo ad un maggior numero di 
particelle espulse dal nucleo con elevata energia cinetica. Segue da cid che 
gli spettri calcolati secondo la relazione (25) includono anche una parte delle 
tracce generalmente considerate come grige e precisamente quelle dovute a 
particeile il cui potere ionizzante specifico si aggira intorno a 4-5 volte il 
minimo. 

Nella fig. 8 sono riportati i dati sperimentali di G. BERNARDINI (!) e L. 
GERMAIN (2°) ottenuti con protoni di 350-400 MeV su emulsioni [ford G5 

Valutando con BERNARDINI l’effetto della trasparenza nucleare e tenendo 
conto che i nostri spettri per energie di eccitazione del nucleo inferiori a 
180-200 MeV si estendono sino a circa 30-35 MeV, abbiamo considerato sola- 
mente quelle tracce di BERNARDINI che si riferiscono a particelle la cui energia 
cinetica non supera questo limite. Ci & stato cosi possibile risalire dai dati 


G. BERNARDINT, E. T. Boorx e S. J. LINDENBAUM: Phys. Rev., 80 905 (1950). 
L. Germain: Phys. Rev., 82, 596 (1951). 

L. Baumnorr e W. Hornine: Phys. Rev., 75, 370 (1948). 

G. BERNARDINI, E. T. Bootu, 8. J. LINDENBAUM: Phys. Rev., 88, 669 (1951). 
J. Fuyimoro e J. Yamacucur: Prog. Theor. Phys., 5, 787 (1950). 
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sperimentali di questo Autore al numero medio di particelle che possono rite- 
nersi, secondo la presente descrizione, di natura evaporativa e quindi risalire 
all’energia di eceitazione termica del nucleo. Risolvendo l’equazione (25), dopo 


30 40 


E,= 250 Mev £,= 450 MeV 


30 
20 


20 


10 


E,= 650 MeV 


© 


Fig. 7 a, b, c. — Spettri energetici delle 
particelle evaporate di carica unitaria 
per differenti valori dell’ energia di 

eccitazione iniziale B,. 10 


aver posto il valore trovato N, = 2,8, si ottiene che in stelle prodotte da pro- 
toni di 350-400 MeV l’energia spesa nel solo processo evaporativo si aggira 
intorno ai 110-120 MeV. 

Per accertare la legittimitä del valor medio N, — 2,8 abbiamo calcolato 
sulla base di esso il numero di stelle in funzione del numero di rami suppo- 
nendo valida la legge: 


— Be = —k il 
(29) PEN a m) = PF, NI, 
1 a) 


’ 

in cui P(N,,N) & la formula di Poisson. I] fattore correttivo LE tiene conto 
del numero limitato m di osservazioni e tende ad aumentare il numero di stelle 
con pochi rami ed a ridurre quelle con molti. Risultati analoghi si ottengono 
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con i dati di L. GERMAIN supponendo che la distribuzione energetica delle 
particelle emesse sia la stessa di quella data da BERNARDINI. 


Il valore trovato Z,—1102120 MeV, 
pur tenendo conto delle notevoli 
estrapolazioni nella stima del numero 
di particelle incluse nell’intervallo 
energetico considerato, non sembra 
essere del tutto in contrasto con le 
recenti previsioni teoriche di BERNAR- 
DINI e coll. (2), qualora si osservi che 
esse sono basate sull’ipotesi che le 
sezioni d’urto differenziali nucleone- 
nucleone nel centro di massa siano 
costanti e che sia lecito trascurare il 
coefficiente di riflessione della barriera 
nucleare. Inoltre & pensabile che la 
stessa adozione del modello di Gold- 
berger favorisca troppo, in quanto 
trascura le transizioni tra stati occu- 
pati e non occupati prodotte dalle 
interazioni nucleari fra i nucleoni che 
costituiscono il nucleo a gas dege- 
nere di Fermi, quei processi che com- 
portano una espulsione di particelle 


20 


stelle in & 


10 


A 


numero di rami À 


0 2 4 6 8 


A G. BERNARDINI e coll. 
A L. GERMAIN 
— Curva teorica 


Fig. 8. — Distribuzione teorica (N; = 2,8) 
e dati sperimentali di G. BERNARDINI € 
L. GERMAIN ottenuti con protoni di 


- : 5 : b 350-400 MeV. 
per urto diretto dei nucleoni primari 


e secondari rispetto a quelli che coin- 
volgono piceoli trasferimenti di energia e che danno luogo al processo 
evapcorativo. : 

In conclusione, il quadro prospettato sembra nelle sue linee essenziali abba- 
stanza soddisfacente, in quanto non solo interpreta in modo corretto i dati 
sperimentali del gruppo di Bristol, ma permette anche, poichè la (25) prevede 
un numero di rami ionizzanti di natura evaporativa minore di quello ottenuto 
da altri Autori (7) (2!) (24), di considerare come coesistenti i due processi: 
a) collisione diretta o indiretta, b) eccitazione ed evaporazione che danno ori- 
gine alle stelle. 


Ci & gradito ringraziare i proff. N. DALLAPORTA, E. CLEMENTEL € G. Pupri 
per le utili discussioni e per l’interessamento con cui hanno seguito questo 
lavoro e la dott. S. CHERSI che ei ha coadiuvato nell’esecuzione dei calcoli 
numerici. 


(4) D. H. Perkins: Phil. Mag., 41, 138 (1950). 
(25) D. H. PERKINS, J. B. HARDING e S. LATTIMORE: Proc. Roy. Soc., 196, 325 (1949). 
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SUMMARY 


The Fermi gas model has been extended to high nuclear temperatures according 
to a mathematical procedure independent of the restrictive condition r/p<1. The 
modification of the momentum distribution brought about by the nuclear interactions 
(Fig. 1) has been considered qualitatively of the same kind as that due to thermal 
agitation. Calculations have been carried out on the basis of the «nuclear equivalent 
temperature » 7 at which the nucleon gas, having a true temperature T, may be con- 
sidered as an ideal one. The relation between true and equivalent temperature has 
been assumed of the form 


(6) t= T + ay(r), 


where a = 0,29 is the Watanabe parameter and y(r) is the Fermi function character- 
istic for higly degenerated gases. The total excitation energy versus true nuclear 
temperature 7 is represented in Fig. 2. On the basis of this curve and following 
Weisskopf’s theory, the evaporation process, initiated by a nucleous of energy Ey, 
has been studied stage by stage taking into account the diminution of the nuclear 
temperature and the variation of the atomic mass number. It has been found that 
the nuclear energy losses can be well represented by a third order difference relation. 
The energy summation formula thus derived leads to a fourth degree equation, the 
solution of which gives the total number of particles evaporated from a nucleous of 
initial energy E, (Eq. (22)). On evaluating the ratio ®(H,) between the number of 
charged particles and the total number N(E,) of them as a function of the excitation 
energy of the residual nucleous, one obtains from (Eq. (22)) the number of ionizing 
particles evaporated from a nucleous of energy H,: 


(25) q(Ey)N? + PEN: — Bo = 0, 


where the functions g(#,) e p(2,) (Eq. (26)) are represented in Fig. 5. Since the mean 
kinetic energy evaluated by means of the energy-temperature curves is sensibly larger 
than that obtained by other Authors, the formula (25) includes black and also some 
grey prongs, as it is shown by the distribution curves (Fig. 7 a, b,c). Therefore, it may 
be compared with the Bristol experimental relation gN?+ pN,— E, = 0, where q = 4 
and p = 37; it is seen that the theoretical previsions are well fitted with experimental 
data (Fig. 6). Fluctuations have been calculated with the Poisson formula. The curve 
represented in Fig. 8, which is in satisfactory accordance with Bernardini’s experi- 
mental evidence, corresponds to N, = 2,8. It is therefore possible to determine by 
the aid of Eq. (25) the energy lost in the evaporation process, which is found to be 
110-120 MeV. This value seems to be not entirely contradictory with theoretical pre- 
visions of Bernardini especially if one considers that the Goldberger model assumed 
by him is based on an unperturbed Fermi distribution. Furthermore, it is believed 
that the evaporation processes are somewhat underestimated by this Author because 
he considers as constant the nucleon-nucleon differential cross-section in the e.m. and 
neglects the reflexion coefficient of the nuclear barrier. 
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Riassunto. — E stato determinato lo spettro differenziale in range dei 
mesoni u, al livello del mare fra 1 e 3,5 GeV/c allo scopo di mettere in evi- 
denza wna supposta anomalia di questo. Dai risultati ottenuti si conclude 
che in questo intervallo di energie non vi sono irregolaritä superiori a 
qualche percento. 


1. — Esperienze effettuate in camera di Wilson da BLACKETT (!) e da 
GLASER (2) hanno indicato la possibilità che sussista una anomalia dello spettro 
differenziala «in momento » delle particelle eariche, al livello del mare. Questi 
Autori hanno trovato infatti che tale spettro differenziale presenta un secondo 
massimo tra 2 e 4 GeV/c con uno scarto di circa il 20% rispetto all’andamento 
medio dedueibile dagli spettri integrali «in range » dei mesoni u. Per quanto 
Vanomalia riscontrata non sia molto significativa a causa degli errori stati- 
stici di quelle misure, ci siamo tuttavia proposti di stabilire se essa si ritrova 
anche in uno spettro differenziale in range. A questo scopo abbiamo realiz- 
zato un dispositivo che permette di determinare, per mezzo di assorbitori, sia 
lo spettro integrale che quello differenziale; questo ultimo mediante misure 
su bande energetiche anche più strette di quelle scelte dagli Autori citati. 
Le misure sono state estese fino a 3,5 GeV/e e sono durate un tempo sufficiente 
ad assicurare, per ogni punto sperimentale, un errore non superiore al 3%. 


G) P. M. 8. Braokert: Proc. Roy. Soc. A 159, 1 (1937). 
() D. A. GLASER, B. HAMERMESH e G. SAFONOV: Phys. Rev., 80, 625 (1950). 
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Operando con questi criteri pensiamo di esserci posti in condizioni favorevoli 
per mettere in evidenza una eventuale anomalia dello spettro, anche se poco 
rilevante o contenuta in uno stretto intervallo energetico. 


2. — Il dispositivo sperimentale da noi usato si compone di due sistemi 
identici, spostati l’uno rispetto all’altro, verticalmente, di 20 em, ciascuno 
dei quali da misure integrali e differenziali. L’insieme & schematizzato in fig. 1; 
i piatti di contatori abex costituiscono il 
primo sistema e i piatti bedß il secondo. 
Tra i contatori e al disopra di essi sono 
posti gli assorbitori di ferro 8, S’ ed A; 
S ed S’ hanno lo spessore fisso di cirea 
20 cm, mentre lo spessore di A & stato 
variato durante le misure, tra 80 e 
300 cm. La frequenza degli eventi regi- 
strati dalle coincidenze tra i contatori 
a,b ec, per i diversi spessori di A, ci da 
la misura integrale; le stesse coincidenze 
abe che non sono anticoincise dai con- 
tatori & rappresentano il numero di 
particelle con un range compreso tra A 
e A+ 8 e danno quindi la misura 
differenziale. La larghezza della banda 
energetica à definita dallo spessore del- 
l’assorbitore S e compreode tutti i me- 
soni con un momento residuo inferiore 
a 300 MeY/c. Cio corrisponde ad una 
banda, all’ingresso nell’assorbitore A, di circa 250 MeV/e per la zona dello 
spettro studiata. Analoghe misure integrali e differenziali sono date dal se- 
condo sistema di contatori nel quale l’assorbitore che stabilisce il taglio per 
la misura integrale à rappreseatato da A + 8, mentre la misura differenziale 
è data dai mesoni che muoiono tra À + 8 ed À + 8+ 8’. I contatori dei 
piatti a, b, ¢ e d hanno le dimensioni di em 4 x 20, mentre quelli delle anticoinci- 
denze « e ß sono di cm 5x40; con la disposizione geometrica da noi adottata 
si ha per entrambi i telescopi un angolo solido di 2 2,3-10-? steradianti. Il 
potere risolutivo del circuito di coincidenza era di 10°5s e l’efficienza delle 
anticoincidenze del 99,8%. Per ragioni di stabilità e per una migliore scher- 
matura dagli sciami laterali, il dispositivo era collocato in un pozzo apposi- 
tamente costruito nei sotterranei dell’Istituto di Fisica di Bologna. 


Fig. 1. — Dispositivo sperimentale. 


3. — Le misure sono state fatte aumentando lo spessore dell’assorbitore A 
da 80 em fino al massimo di 300 cm e diminuendolo successivamente fino a 
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ritornare al valore iniziale, Il buon accordo tra i valori ottenuti in queste due 
fasi della misura da entrambi i dispositivi ci ha permesso di concludere la non 
esistenza di derive sistematiche strumentali. 

Per avere un controllo continuo del regolare funzionamento dei due dispo- 
sitivi lo spessore di A é stato variato sempre di 10 o 20 cm per volta; tenendo 
presente lo spessore di S e di 8’ si vede che, in tal modo, i due dispositivi 
danno, sia per lo spettro integrale che per quello differenziale, delle serie di 
punti sperimentali direttamente confrontabili; l’accordo tra i valori ottenuti 
dai due dispositivi é risultato sempre ottimo e le misure sono state quindi 
sommate migliorandone la precisione statistica. 

Abbiamo voluto anche stabilire l’influenza che avrebbe potuto avere lo 
scattering multiplo nell’assorbitore sulla sensibilità della misura differenziale. 
Il calcolo della distribuzione gaussiana dell’angolo di scattering multiplo all’en- 
trata dello strato S o 8’ di fine range, nell’ipotesi di v/e 21, si pud effet- 
tuare (?) con la: 


Tr 


a ‘de ER [de Fifty 2 1 | 
1X, | EX, | (Eye) "X, \Hi—ym E—ym 3 


ey 


TABELLA I. TABELLA II. 

| ae | GeV/e | Media oraria | | Sem: | Siem: GeV/e | Media oraria 
Fe arıa | Be aria | | 

ee 2 # 2 Kl SET 
200 21320.160,72, 16.173.072 LT 585 ATO Lie it 13,3 + 0,6 

| 520 | 420 | 0,92 | 155,.0+16 || 715 5701) 11823. | 1122208 

650 | 590°) CTI \ 143,62 10 || 845 | 670 | 1,45 | 10,5 +02 
780 620 | 1,35 | 130,9 + 0,84 975 770 | 1,67 | 9,66 + 0,2 
910 | 720 | 1,57 | 121,6 + 0,94 | 1105 870 | 1,90 | 8,27 +0,17 | 

| 1040 2204 MTS) RIS 7 OS | 1.2235 | 9700! 2,18 7,61 + 0,17 
1 170 920 | 2,05 |-103,8+0,7 | 1365 | 1080 | 2,45 | 6,93 +0,14 | 
1300 | 1030 |.-2,30 96,5 + 0,6 1495 | 1180 | 2,68 6,23 + 0,11 
1430 | 1130 | 2,55 | 89,2+0,6 | | 1560 | 1230 | 2,80 | 5,98 + 0,14 
1295 HET 1607) ee, TV 8898.08 1625 | 1280 | 2,93 | 5,68 + 0,11 
1560: | 1230 | 2,8 | 34.0408 | 1690 | 1330 | 3,03 5,29 + 0,14 
1625 | 1280 | 2,9 82)8'4-.0;7, | 1755 | 1380 | 3,08 | 5,10 + 0,09 

| 1693 | 1330 | 3 79,4 + 0,4 | 1810 | 1420 | 3,12 4,83 + 0,11 

175° 1,370°| 3,10 | 27,5 + 0.5 1862 | 1460 | 3,17 | 

| 1797 | 1410 | 3,21°|, 75,240,4 | | 1914 | 1500 | 3,42 4,75 + 0,14 | 
1849 | 1450 | 3,30 | 72,8 + 0,6 1965 | 1540 | 3,54 4,45 + 0,14 
1901 | 1490 | 3,40 70,8 + 0,6 


(3) B. Rossı e K. GREISEN: Rev. Mod. Phys., 13, 264 (1941). 
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Run 


dove Hy — 21 MeV; X, — 14,4 g/cm? é il cammino di radiazione nel ferro; 
E,— yx, ed E,— yx, sono l’energia rispettivamente all’entrata e all’uscita del- 
l’assorbitore A, nell’ipotesi che la perdita y di energia per g/cm? sia costante. 

Per mesoni di energia 0,9 GeV e 3,2 GeV (estremi della parte dello spettro 
da noi esplorata) e per una energia finale di 0,22 GeV, che rappresenta la lar- 
ghezza a fine range della banda energetica da noi scelta per la misura differen- 
ziale, si ottiene: 


65.2159 e 6,,213°. 


Considerando il valore di questi angoli si conclude che la diffusione della 
banda, per effetto dello scattering multiplo, si mantiene di entità modesta 
rispetto alla larghezza della banda stessa e percid lo scattering non pregiudica 
la sensibilita della misura agli effetti della rivelazione di una anomalia che 

dovrebbe occupare una regione 
108g" ls"! sterad 107? cm’? s"\.sterag”! ie dello spettro che é circa quat- 


9 


to tro volte la larghezza della 
à Ie banda medesima. 

F Le tabelle I e II riportano 
5 a le frequenze rispettivamente 


della misura integrale e di 
quella differenziale. 
5 Ike Per il passaggio dai g/cm? 
di ferro impiegato nelle varie 
misure a quelli equivalenti di 
19% aria & stato usato il nomo- 
gramma di Wick (*) e per il 
successivo passaggio ai momenti 
cy, . ei Siamo serviti delle curve di 
: Princeton (5). I risultati ripor- 
tati dalle tabelle sono gia cor- 
retti dal fondo di casuali e di 
m Spettro integrale di range. sciami e dal contributo, alle mi- 
Le linee a tratto continuo sono le curve di Rossi. sure differenziali, della ineffi- 


g.cm 


Fig. 2. - e Spettro differenziale di range. 


cienza delle anticoincidenze; 

queste correzioni sono state valutate sperimentalmente per vari punti del- 
l’intervallo energetico studiato. 

Altra correzione da apportare alle misure avrebbe potuto essere quella 

relativa all’effetto barometrico. Il coefficiente barometrico, assegnato per la 

radiazione totale incidente, diminuisce certamente al crescere del taglio sullo 


(4) G. C. Wick: Nuovo Cimento, 1, 312 (1943). 
(5) J. X. Montgomery: Cosmic Ray Physics (Princeton 1949), p. 354. 
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spettro; calcolare l’effetto barometrico assumendo per tutto lo spettro un coef- 
ficiente costante uguale a quello di DUPERIER (*) porta quindi a sopravalu- 
tare la correzione per l’intervallo energetico nel quale si sospetta l’anomalia. 
Per le variazioni di pressione atmosferica registrate nel corso delle nostre mi- 
sure, le correzioni, pur assumendo un tale coefficiente costante, erano ampia- 
mente contenute negli errori statistici e percid non abbiamo creduto opportuno, 
ai fini di questo lavoro, di tenerne conto. 

La fig. 2 mostra l’andamento che hanno i nostri punti sperimentali norma- 
lizzati sulle curve di Rossi. 

4. — Dopo quanto é stato esposto nei precedenti numeri & facile concludere 
che una anomalia nello spettro differenziale dell’ordine del 10% non avrebbe 
potuto sfuggire. A conclusione di questo lavoro possiamo quindi dire che non 
esistono anomalie nello spettro differenziale in range dei mesoni u. al livello 
del mare nell’intervallo di momenti da 1 a 3,5 GeV/e di entita superiore a 
qualche per cento. Resta sempre aperta la possibilità che esistano irregolarita 
più modeste che non riguardano perd i fini della presente ricerca. 

Vi & tuttavia un punto che bisogna chiarire: la possibile anomalia citata 
all’inizio, supposto pure che abbia consistenza nei riguardi della precisione 
statistica delle misure, è emersa in uno spettro in momento della radiazione 
carica al livello del mare in camera di Wilson, in campo magnetico; il presente 
& invece uno spettro in range ottenuto per assorbimento, del quale possiamo 
con certezza affermare che esso riguarda esclusivamente mesoni w che sono 
le uniche particelle le quali, a causa della debolissima interazione nucleare, 
sopravvivono sotto notevoli spessori di assorbitori. Anche nello spettro in 
momento, determinato in campo magnetico, & abbastanza plausibile, per quanto 
si sa attualmente sulla radiazione cosmica nell’atmosfera (7), che siano pre- 
senti essenzialmente solo mesoni u, giacchè le componenti assorbite rapida- 
mente per interazione nucleare possono ammontare globalmente soltanto a 
qualche percento del totale delle particelle penetranti. Pur non potendo esclu- 
dere un intervento di queste componenti, sembra piuttosto improbabile di 
potere imputare ad esse l’anomalia in discorso, che resterebbe quindi, in base 
alla nostra esperienza, da escludere. 


5. — Questa ricerca & stata fatta nell’Istituto di Fisica dell’Universita di 
Bologna, al Direttore del quale, prof. G. VALLE, porgiamo i nostri vivi ringra- 
ziamenti per i mezzi largamente messi a disposizione. Egualmente vivi ringra- 
ziamenti porgiamo al prof. A. RoSTAGNI, Direttore dell’Istituto di Fisica del- 


(8) L. JANossy: Cosmic Rays (Oxford, 1948), p. 195. 
(7) G. Puppr e N. DALLAPORTA: Progress in cosmic rays physics (Amsterdam, 1952), 
p. 315. 
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l’Universitä di Padova, che in tutti i modi ha appoggiato e facilitato la nostra 
collaborazione. Al prof. G. Pupp, che ci ha suggerito l’esperienza e ci & stato 
largo di preziosi incitamenti e consigli durante la sua esecuzione, rivolgiamo 
l’espressione della nostra più viva riconoscenza. Porgiamo infine i nostri sinceri 
ringraziamenti all’Amministrazione Comunale di Bologna per l’aiuto presta- 
toci, all’ing. L’ABBATE della SABIEM e all’ing. BARCELLONI Corti della TIMO 
per il materiale messo a disposizione. 


SUMMARY (*) 


The range differential spectrum of mesons u has been determined at sea level 
between 1 and 3,5 GeV/c in order to gain evidence of an anomaly it is supposed to 
present. From the results obtained Authors conclude that in this energy interval 
there are no irregularities greater than some percent. 


(*) Editor’s translation. 
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Misura della concentrazione in deuterio dell’acqua di condensa 
di un soffione boracifero di Larderello. 


E. Cerrar, ©. MARCHETTI e M. SILVESTRI 
Laboratori CISH - Milano 


(ricevuto il 21 Aprile 1952) 


Riassunto. — Mediante il metodo del galleggiante e con un prearricchi- 
mento noto, & stata confrontata la concentrazione in D,O dell’acqua di 
condensa di un soffione di Larderello, con quella dell’acqua di rubinetto 
di Milano. Le due concentrazioni risultano eguali, entro qualche percento. 


1. — La concentrazione in deuterio dell’acqua naturale (espressa in con- 
centrazione molare di D,O) & di circa 0,015%, variabile leggermente a seconda 
della provenienza dell’acqua (1). Abbiamo assunto che l’acqua di rubinetto di 
Milano abbia concentrazione in D,O pari a quella indicata e ci siamo proposti 
di confrontare con essa la concentrazione di D,O dell’acqua di condensa in 
un soffione boracifero della zona di Larderello (Toscana). 

Per la misura della concentrazione in deuterio dei campioni siamo ricorsi 
ad un metodo densimetrico molto noto, quello del galleggiante, nella variante 
suggerita da SAPIRSTEIN (?), utilizzato come metodo assoluto (?). Questo con- 
siste nella determinazione della temperatura f, a cui un galleggiante di vetro 
di opportune dimensioni si stacca dal fondo di una piecola ampolla contenente 
acqua naturale distillata, o comunque ad una concentrazione isotopica nota, 
e nella successiva determinazione della temperatura it a cui si distacca lo stesso 
galleggiante quando l’ampolla contiene l’acqua a concentrazione ignota. Questa 


(1) I. KIRSCHENBAUM: Physical properties and analysis of heavy water (New York, 
1951), p. 398. 

(2) L. A. SAPIRSTEIN: Journ. Lab. Olin. Med., 35, 793 (1950). 

(8) E. Cerrar, C. MARcCHETTI e M. SıLvEsırı: Nuovo Cimento, 9, 530 (1952). 
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concentrazione per un dato galleggiante (cioè per un dato t,) & funzione della 
grandezza t—t,(Atı,). La sensibilità di questo metodo con l’apparecchiatura 
da noi usata si aggira intorno ad 1/20000 in concentrazione molare di D,O 
nell’H,O, operando con volumi di liquido pari a circa 2 cm. 


2. — Di solito ogni campione di acqua, sia naturale che arricchita, viene 
sottoposto ad una purificazione che consiste in una néutralizzazione con so- 
stanze non contenenti acqua o idrogeno, quando occorre, ed inoltre ad un 
trattamento con pochi mg di KMnO,. Questo è seguito da tre distillazioni 
sotto vuoto, dopo di che si pud procedere alla misura usando un piccolo bagno 
termostatico corredato da un termometro al decimo di grado, letto mediante 
un cannocchiale. In questo caso perd, dovendosi trattare un’acqua che con- 
tiene fra l’altro CO,, H,S ed NH, diseiolti, si e proceduto ad una puri- 
ficazione particolare. La CO, non costituisce una preoccupazione in quanto 
per semplice degasaggio é facilmente eliminabile, mentre gli altri due composti, 
se non bene eliminati, in seguito ad ossidazione dell’H,S avrebbero potuto 
fornire piccole quantita di (NH,),SO, a scapito della purezza del campione, 
e quindi della precisione della misura. Tutta l’acqua di condensa da noi usata 
per la misura é stata percid trattata nel seguente modo: addizionata con car- 
bonato, anidro e secco, di sodio e potassio e permanganato potassico, pure 
secco, per spostare la maggior parte dell’ammoniaca ed ossidare i solfuri, é 
stata bollita all’aria per aleuni minuti. Quindi per neutralizzare l’alcalinità e 
l’ammoniaca residua é stata addizionata con pochi cristalli di pirosolfato sodico 
secco fino a reazione leggermente acida. Poi per eliminare i sali disciolti il tutto 
é stato accuratamente distillato a pressione normale. Un campione di acqua 
cosi ottenuto & stato sottoposto alle solite tre distillazioni sotto vuoto e quindi 
alla misura col metodo prima descritto. 


3. — Poiche, come abbiamo detto, la precisione del metodo è di circa 
1/20000 in concentrazione molare, e noi dovevamo misurare dei valori nel- 
l’intorno di 1,5/10000 (concentrazione naturale) abbiamo eseguito una misura 
preliminare usando anziché un termometro al decimo di grado, un termo- 
metro di Beckmann al centesimo in modo che fosse apprezzabile qualche mille- 
simo di grado. 

La misura preliminare ha dato il seguente risultato: 


a) Acqua naturale. — Lettura per tre volte di seguito sul termometro 
Beckmann della temperatura t,, a cui il galleggiante si distacca dal fondo del- 
l’ampolla che contiene il campione di acqua naturale: 


to = tp +0,460°C; tt — tp + 0,455 00; t — 13 + 0,455 0; 
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mediando si ha: 


ibe == tp a 0,457 0C . 


b) Acqua Larderello. — Lettura per due volte di seguito sul termometro 
Beckmann della temperatura ¢ a cui il galleggiante si distacca dal fondo del- 
Vampolla che contiene il campione di acqua in esame: 


t = ty + 0,457 °C; t = t, + 0,459 00; 
mediando si ha: 
t = te + 0,458 °C. 
Mt == tty — 0,001.00": 


Questa misura non rivela alcuna differenza apprezzabile fra le concentra- 
zioni isotopiche delle due acque. Il metodo, infatti, per ogni differenza nelle 
concentrazioni pari all’1% fornisce un valore di At ch’è intorno a 3 0 3,5 °C. 
In tal modo, tenendo conto di tutte le possibili cause di errore, di fronte alle 
quali la precisione di lettura della temperatura mediante il Beckmann diviene 
sovrabbondante, si pud raggiungere una sensibilita di poco superiore ai valori 
stessi delle concentrazioni che vengono misurate. Per questa ragione siamo 
ricorsi ad un metodo che amplificasse le eventuali differenze non rivelate da 
questa misura. I] metodo consiste semplicemente nel far subire un certo arric- 
chimento isotopico ai campioni da confrontare, facendo loro subire una deter- 
minata elettrolisi. 


4. — L’arricchimento dei campioni viene effettuato in una celletta verticale 
in acciaio inossidabile con elettrodi cilindrici coassiali, di cui quello interno 
& raffreddato ad acqua corrente (fig. 1). Si usano soluzioni alcalinizzate con 
quantitä note di KOH e si pesa il quantitativo iniziale e finale di ogni elettro- 
lisi. Si opera su quantitä pari a circa 100 g che vengono ridotti di 15 o 20 volte 
a corrente costante ed in modo che sia minimo lo schiumeggiamento e nullo il 
trasporto di elettrolito nei gas. 

Si sono eseguiti in tempi successivi, nell’ordine i seguenti arricchimenti: 


a) g 100,04 di acqua naturale con aggiunta di g 1,39 di KOH sono stati 
elettrolizzati fino a g 7,5 di acqua arricchita. 


b) g 100,00 di acqua Larderello con aggiunta di g 1,39 di KOH sono 
stati elettrolizzati fino a g 5,72 di acqua arricchita. 


c) g 100,00 di acqua naturale con aggiunta di 1,39 g di KOH sono stati 
elettrolizzati fino a g 6,40 di acqua arricchita. 
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Il tutto & avvenuto con una corrente costante di 4 + 0,2 A ed alla tem- 
peratura di 16 + 0,5 °C. 


Fig. 1. 
1. Tubo ingresso acqua raffreddata 
(Mat. ottone) & 6x3. 


2. Tubo uscita acqua raffreddata 
(Mat. ottone) @ 6x3. 


3. Coperchio in plexiglas 
(Mat. plexiglas) @ 58, spess. 15 mm. 


4. Elettrodo esterno ausiliario 
(Mat. ace. Inox) @ 49x45. 


5. Tondino di sospensione dell’insieme 
(Mat. acciaio trafilato) filo @ 1,5. 


6. Fascia per la sospensione dell’insieme 
(Mat. acciaio dolce) & 60x50 


7. Elettrodo intermedio 
(Mat. acciaio Inox) & 30x28. 


8. Elettrodo interno con raffreddamento 
(Mat. acciaio Inox) & 22x20. 


9. Distanziatore in politene 
(Mat. politene). 


I campioni raccolti, neutralizzati con AgNO, secco e distillati sotto vuoto, 
hanno dato i seguenti risultati (si è usato il termometro al decimo letto col 


cannocchiale e si sono eseguite tutte le misure in una unica mandata): 


TER. i lee bs 

a) Acqua naturale arricchita: à = 30,10 | t = 30,105 
ty = 29,80 & 

b) Acqua naturale normale: fp = 29,805 | to — 29,803 
Sy) t = 30,18 1: 

c) Acqua Larderello arricchita: — 30,175 7 290.178 
ie Sl oa me 

d) Acqua naturale arricchita: — 30,14 | t — 30,145 
t — 29,80 | ke 

e) Acqua naturale normale: fy = 29,80 | to = 29,800 
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Come acqua naturale normale di riferimento, ei si è riferiti alla misura più 
vicina a quella eseguita sull’acqua arricchita. si & ottenuto percid: 


a’) Acqua naturale arricchita: Atto = 30,105 — 29,803 = 0,305,,, e035 
la concentrazione molare corrispon- 
dente è: 7% = 0,104. 

b') Acqua Larderello arricchita: Att) = 30,178 — 29,803 = 0,3752¢, s033 
2% = 0,124. 


e') Acqua naturale arricchita: Atto = 30,145 — 29,800 = 0,345, 800: 
29%, 0,116. 


I due confronti con l’acqua naturale trattata prima e dopo dell’acqua di 
Larderello sono stati eseguiti per attestare la riprodubciilita del metodo entro 
il limite previsto degli errori. 

Poiché sono stati eseguiti arricchimenti continui, a corrente ed a tempe- 
ratura costanti e ©, ©, <1, si pud applicare la formula: 


V PACE 
ee 


il volume finale di acqua a concentrazione x in deuterio, 


in eui: 


NX 
& 


V, è il volume iniziale di acqua a concentrazione x, in deuterio, e 


q il fattore di separazione della cella elettrolitica. 


Usando in luogo dei volumi i pesi, si commette un errore trascürabile ; 
introducendo poi nella formula scritta nel modo seguente: 


i dati riguardanti l’esperienza a’) sull’acqua naturale, sapendo che ©, = 0,015% 
(concentrazione naturale) e risolvendo rispetto a q, si ottiene: q — 3,85. 

Se si risolve ora la stessa equazione usando le grandezze dell’esperienza c’) 
sempre sull’acqua naturale (e quindi ancora con % = 0,015%) si ottiene: 
q = 3,74. 

I due valori di q cosi ottenuti concordano bene entro eli errori sperimentali 

Abbiamo assunto allora quale valor medio del fattore di separazione della 
cella, il valore Qmedio = 3,8. 

Se ora, noto questo valore, si risolve nuovamente la stessa equazione con 
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i dati della esperienza 5b), eseguita sull’acqua di Larderello, questa volta ri- 
spetto alla concentrazione iniziale, si ottiene: æ, = 0,0149%. 

Stimando che gli errori sperimentali assommino a circa 4% (2% per l’in- 
certezza del valore di q ottenuto entro il 5% e circa 3% per l’errore di lettura 
del termonetro) si pud concludere che il valore della concentrazione molare di 
deuterio sull’acqua di condensa del soffione è: x, — 0,0149 + 0,0006%, cioè 
praticamente risulta uguale, entro meno dell’1%, all’acqua di rubinetto della 
città di Milano. 


SUMMARY 


Heavy water concentration of the Larderello (Toscana) natural water vapor has 
been measured by a float method. In order to rise precision, the sample was previously 
enriched in an electrolitic cell, whose separation factor had been controlled by Milan 
tap water, before and. after the enrichment of the Larderello sample. Assuming the 
concentration of the former as 0,015%, the Larderello water vapor concentration 
results 0,0149 + 0,0006 %. 
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On the Energy Range Relation for fast Muons in Rock. 


M. Manno and L. RONCHI jr. 


Institute of Physics of the University - Arcetri, Firenze 


(ricevuto il 22 Aprile 1952) 


Summary. — A discussion is given of the best estimate of various types 
of energy loss; the energy-range relation for muons in rock and water 
is then evaluated. Apart from the uncertain amount of «nuclear » losses 
this relation should be quantitatively correct to about 1%. Finally two 
different curves are given, corresponding to zero «nuclear» losses and 
1.0 x 10-8 g-! cm? fractional «nuclear » energy losses respectively; within 
them, linear interpolation (for any amount of « nuclear » losses) is allowed. 


1. — Introduction. 


A fairly accurate knowledge of the energy-range relation for muous in 
rock is most important for the interpretation of any cosmic ray experiment 
underground. 

No such relation, among those, which have been used so far, can be relied 
upon quantitatively; the present paper is intended as a step towards a com- 
plete solution of the problem, through a detailed discussion of the various 
difficulties, the improvement of some approximate evaluation previously 
adopted by other authors and a final evaluation of the range by taking into 
account only collision, radiation and pair creation losses, which for brevity’s 
sake will all be indicated henceforward as «extranuclear losses » We believe 
that this evaluation, within its limits, is fairly accurate; we hope that it may 
help to reduce the existing uncertainties in the comparison between theory 
and experiments, and may thus allow to give experimentally stricter limits 
for the amount of other Josses. 

To this purpose the energy-range relation will also be evaluated for an 
assumed «nuclear» loss of 1.0x10-* g-!cm?; linear interpolation (and pos- 
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sibly moderate extrapolation) will thus allow an easy comparison with future 
experiments. 


2. — Collision losses. 


Since we are primarily concerned with very high energies, of the order of 
100 GeV, we can assume, for the range of a 5 GeV kinetic energy muon, an 
average value of, say, 25 hg/cm? (1) and confine ourselves to the energy loss 
of particles with more than 5 GeV kinetic energy. This introduces a very 
small error (?), which, if the case may be, can be easily corrected and allows 
many simplifications to be made, for instance to dispense with the exact 
knowledge of the chemical composition of the rock, as we shall see. 

Since there are some differences among different authors, we reproduce 
here our starting formulae. 

For muon collision probability against electrons we use the formula given 
by BHABHA (*) and Massey and ÜORBEN (§) for 1/2 spin particles, namely: 


2nNZ „me: QE" ee E' \ 
A ‘fr En 


il (BZ, EN AE = —— rn + < lu. — 
() A ) j Eu 2\E+ pe? 


(cfr. also Rosst and GREISEN (§), page 244), where Æ is the energy of the 
muon (4 its rest mass), H’ the energy transferred to the electron and PY the 
maximum transferable energy. 

Formula (1) holds for sufficiently large values of EB’, say E> », for which 
electrons may be considered as free. By multiplying (1) by #’, and integrating 


in the range 7 < E’< EY, we get 
db 2n NZ mel, E, 1 ER 

‘ 2 m do m 

(2) er ee re a ap er 
dx Joon > » A B N 4 (E + ue)® 


allowance being made for the fact that in our case n may be chosen <E,. 


(1) We use hg/em? as a unit, in order to conform with the widespread usage af 
meters w.e.; by our notation, however, we hope to avoid possible confusion by empha- 
sizing the purely massie equivalence. 

(2) For instance, Wick’s value ($) of the range of 5 GeV muon in iron is 28.3 ho/om?, 
in water 21.6 hg/em?; Montgomery’s (D. J. X. MONTGONMERY: Cosmie Ray Physies, 
Princeton 1949, p. 353) value is 22 hg/em® in air. 

(3) G. C. Wick: Nuovo Cimento, 1, 302 (1943). 

(4) H. J. BuABHA: Proc. Roy. Soc., 164, 257 (1938). 

(5) H. S. W. Massey and H. C. CoRBEN: Proe. Cambr. Phil. Soe., 85, 463 (1939), 

(°) B. Rossı and K. GREISEN: Rev. Mod. Phys., 18, 240-309 (1941). 
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In the interval 0 < £’ <7 we will adopt Bethe’s (*) (§) loss formula, na- 
mely: 


2aNZ mel 2me®n 
3 —— = > In ——__ — 57), 
- | Ei: 7 A F p? | (1 — p*)I*(Z) ® 


The total loss is therefore: 


dE 2aNZ _ me 2mv°E! re | Er 
(4) ( ) = à ae À, 
G2 Jeon A BP 1— p2)I*(Z) 4 (E+ ue)?] 


where 1(Z) is the average ionization energy of the atom considered. 
For the density effect correction we adopt Fermi’s (*) limit formula (v— €) 
[formula (39)], which can be written as follows 


(5) A=— B) 


| 


2aNZ me? pm Ne Z 
Ao (1— B2)12(Z) am ~ A 


so that the final corrected expression for collision loss comes out to be: 


ameH. 2A 1 EE 


1 
6 RE; (5 = i | 
(6) Neh? 0Z PTimr uc?) (E + ne)? ' 


| = 2aNZ . me 
coll ” 


de, A "0 & | 


At the high energies we are considering, formula (5) should be accurate 
enough; as a matter of fact, an approximate calculation based on HALPERN 
and Harr’s (2°) method, shows that, in the most unfavourable case (iron) the 
relative error in formula (6) may amount to about 1% at 5 GeV and decreases 
obviously with increasing energy. 

Had we used, instead of (3), a formula given by BLOCE (7!) and adopted 
by HEITLER (#2) and HALPERN and HALL (**), the final energy loss would differ 
by some percent; we adopt (3) following Rossi and GREISEN (‘), FERMI (*), 
and Wick (*), who has checked it (cfr. also a discussion of this point by BETHE 
and LIvINGSTON (§)). Other correction terms, also given by BETHE and Lr 
VINGSTON (8) are altogether negligible at high energies. 

Finally it must be emphasized, as shown by FERMI (*) and, with more 
detail, by SCHÖNBERG (#*), that Cerenkov losses are included in our final for- 


(?) H. A. Berne: Ann. d. Phys. 5. 325 (1930); Zeüs. 7. Phys., 76, 293 (1932). 

(8) M. S. Lrvryeston and H. A. BETBE: Rer. Mod. Phys. 9, 245 (1937). 

() E. Ferm: Phys. Rev. 57. 485 (1940). 

(2) O. Harper and H. Hatt: Phys. Rev. 73. 477 (1948). 

(2) F. Broca#: Zeiis. j. Phys., 81. 363 (1933). 

(2) W. Herrzer: The Quantum Theory of Radiation (Oxford, 1944, reprinted 1949), 
page 218. 

(3) M. SCHÖNBERG: Nuovo Cimento. 8, 159 (1951). 
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mula. Quite recently SCHÖNBERG (14) has cast some doubt on the correctness 
of the classical treatment, on which the calculation of the density effect cor- 
rection is based; but, in any case, he concludes rather for a different distri- 
bution of the total loss between Cerenkov and other collision losses than for 
a different total amount of loss; the experimental evidence, as far as it goes, 
seems rather to confirm the validity of the formulae on which our formula (6) 
is based. Hereupon we will confidently use it. 


8. — Radiation losses. 


For the differential radiation probability by muons we will adopt the same 
formula as CHRISTY and KUSAKA (#5), namely 


16 -__ Z2/m\? (3 1—v 121—v) # 1 
Ane ce y2 phe ME EA 
(7) D(H, v) dv 3 aNrè 7 (2) (= He ; )| D ue 2 


v being the fractional energy loss and # the energy of the muon. When 


multiplied by Hv and integrated, formula (7) gives for the average energy 
loss the expression (1) 


dH Z [m\? E 1g 1 
on = 4«N? —|—) # |In — + In — —-|. 
Se | dx i ae A Bi | | 2 ue? BZus A 


Pa 


which we will use for the calculations in n. 7. 

We may remark that the application of range concept and average energy 
loss evaluation is open to doubt for radiation losses, but we will not discuss 
this point here, because the contribution of radiation losses to the total is 
only a small fraction for most part of the range. 


4. — Pair creation losses. 


We have recently (17) worked out a formula for pair creation losses, namely 


d N s \2 Ds 
7) REIN lan. 
de Jpair Au 7 ue” 


(9) 


(4) M. SCHÖNBERG: Nuovo Cimento, 9, 210 and 372 (1952). 

(5) R. F. Curısty and §. Kusaxa: Phys. Rev., 59, 416 (1941). 

(6) E. P. GEORGE: Progress in Cosmic Ray Physies, Edited by J. G. Wilson, 
Amsterdam, 1952, page 409, gives the same formula, with the constant term within 
brackets taken as — 1/2; this difference, however, is of no practical importance. 

(7) M. MaxDd and L. Roneout jr.: Nuovo Cimento, 9, 105 (1952). 
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This formula, however, neglects the effect of the screening; to take it into 
account accurately would require a long calculation work with cumbersome 
formulae and has not yet been carried out. We may, however, hope that 
an approximate treatment, as already used by HaYAKawa and TOMONAGA ('°) 
and based on the extrapolation of Bhabha’s (#) limit formulae, though not 
giving completely carrect results, may at least supply us a reasonable estimate 
of the reduction introduced by the screening. 

According to BHABHA (:*) the cross section for the production of a pair 
of total energy E’ by a muon of energy E is given by 


(10) o(E, E')dE'= > (xZr,)*f(E, E')(dE’/E'), 


where f(E, E’) has the following approximate values, if screening is neglected 


fx 
(11) f(E, E') =In M, for 2me<E'<Emjn, 
me uE 
\2 9 
(12) f(E, E') = (iF) m2, for Emju<E'<E. 
\uE' ue > 


If, on other hand, screening is taken into account, it comes out that (11) 
is correct only for 2me<E'’<2ame, where a is a constant of the order of 
aZ-3, while for Zam® < E' <Em/u we should use 

Em 


(13) HE,E') =Inaln aE 


Similarly (12) can be used only for 2(E ue} -me/a< E'< E, while for 
Em/u< E' < 2(E/uc*)?me*/a we should use 


(4) HE, E') = = (=) In (a ne) 


To get the total loss, Hayakawa and TomonaGa (38) integrate (W/A)E’ x 
xX o(E, E’) dE’, by using for f(E, E’) the expressions (11), (12), (13), (14) over 
the ranges of E’ above indicated. Although these are only limit formulae 
which do not even overlap with one another, they give nevertheless a fairly 
accurate integrated value, as we have shown in the above quoted paper (*’). 
Much more so we may expect that such formulae will give a good estimate 


(48) S. Hayakawa and S. TomonaGa: Progr. of Theor. Phys., 4, 287 (1949). 
(1) H. J. BHasBHA: Proc. Roc. Soe.. A 152. 559 (1935); Proc. Cambrigde Phil. Soc., 
31, 394 (1935). 
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for the ratio of the losses as calculated with screening included and without 
screening; by multiplying this ratio by our more accurate evaluation of 
losses (9), we hope to get a quite reasonable value for pair creation loss. 

The losses, as evaluated{by the approximate method of HAYAKAWA and 
TOMONAGA previously outlined, would be: 


N 16 m. 2E| 14 
(«Zr,)? ai | E 5 Ina+1 3 In = 9 we? (1 —In 2)In 


+ 


auc? 


i 


this is, we believe, the correct form of Hayakawa and Tomonaga’s formula 
(also quoted by GEORGE, I. c¢. (1), page 410) which however contains some 
misprints. 

If screening could be neglected we would obtain, by the same method, 
a loss 


+a (yg — 2022+ int 2)—(1—In 2) In = — (221m 2 + Int 2) 


28 


8 N 16. E 
(16) > rar) mir | ne 
IT 


14 m 
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mette) (2— 2 In 2 + In?2) | 
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The ratio may be simplified, by keeping only the mäin terms, as follows 


af Ina+1 
screening oe 
no screening 16, E 14 } 
Rene ae 2 sen 
RC 9 


g 


By combining (9) and (17) we get what we consider as the most accurate 
evaluation to date of pair creation losses, namely 


E 16 
19.: 00 NI | Se 
| A ne nn (199 53.7) (in a+4) 
0e Aus eve 16. EF old * 
— In — — — +1n2 
9 N ue 9 is 


5. — Other losses. 


Until recently it was generally believed that «extra nuclear » losses practi- 
cally exhausted all possible causes of absorption of fast muons. 
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Observations in nuclear plates underground by GEORGE and Evans (2), 
Aowever, have shown that muons may also give rise to nuclear stars; besides, 
@enetrating pairs have also been observed underground by BRADDICK and 
HENSBY (**), Nass (*?) and WOLFENDALE (35); these pairs may also be reaso- 
mably ascribed to some hitherto unknown interaction of fast muons (74). The 
anergy losses corresponding to these phenomena and which, for brevity’s sake, 
will be called henceforward « nuclear » losses, have been evaluated by various 
authors. with widely different results. 

GEORGE (>) estimates that star production is responsible for a fractional 
energy loss of about 6.5 x10-* g cm?, while the fractional loss due to pene- 
trating pair production is tentatively estimated by him as 2.0 X10-5 g-! em. 
GARELLI and WATAGHIN (2%) evaluate theoretically this loss as 3.5 x10-* g-? em?. 
Hayakawa (27), on the other hand, believes that «nuclear » losses have been 
overestimated by the above mentioned authors. 

A conclusive comparison with experiments cannot yet be done, in our 
opinion, because to date both the intensity-depth relation (deep underground) 
mand the sea level spectrum of high energy muons are far from being well 
established experimentally. 

If, anyway. we assume the intensity-depth relation, which has been ela- 
Morated by GEORGE (1. e. (1%), page 415, fig. 9) after a careful analysis of the 
existing data, to be exact (?5), a remarkable agreement is found between 
George’s «average » experimental curve and the theoretical one, also evaluated 


(®) E. P. GEORGE and J. Evans: Proc. Phys. Soe., A 63, 1248 (1950); A 64, 
-193 (1951). 

(@) A. J. J. Brappick and G. Henssr: Nafere. 144 1012 (1939). 

(2) F. Nasm: Ph. D. Thesis Manchester, as quoted in the article by as 
(Ql. e. (5), page 411 and 448). 

ES) A. WOLFENDALE, as quoted by GEORGE, 1. ce. (2%), page 411 and 499. 

(*) Accurate observations at great depths by Coccoxı and others, however, set 
“an upper limit to the cross section for penetrating pair production, which is conside- 
rably lower than the value 4-10-% cm? per nucleon estimated by Nasx amd WorL- 
FENDALE. 

Added in proofs. — The latest results of BRappicK. Nasa and WOLFENDALE (Phil. 
Mag., 42, 1277, (1951) seem to confirm, at least for lead, their high value of the 
-eross-section, obtained with a controlled cloud chamber technique, but AMALDI, Ca- 
"STAGNOLI, GIGLI e Scruri, with conuters and hodoscope technique. find again a much 
lower value. 

G5) E. P. GEORGE, 1. ec. (2%), page 411. 

ES) ©. M. GarerLı and G. WaracHın: Nuovo Cimento. 8. 695 (1950); Phys. Rev., 
79. TIS (1950). 

(7) S. Hayarawa: Phys. Rev. 84 37 (1951). 
ë ES) We may remark that BarNOTHY and Forro’s result at 2550 hg/em? which 
appears to depart most from the average curve, should not have been included because 


ene 2 
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by GEORGE, on the basis of his energy losses and a muon sea level differential 
spectrum of the familiar power law type (43). This agreement, however, 
cannot be taken as a confirmation of the high value (2.5 x10-*#) of «nuclear 
losses » used by GEORGE, unless the exponent in the power law spectrum is 
confirmed. 

Now, as remarked by GEORGE himself, such spectrum is the most natural 
extrapolation from CARO, PARRY and RATHGEBER’s (29) results, provided that 
the sea level spectrum reflects, apart from muon decay and absorption, the 
generation spectrum of pions; but this is only true if we accept the hypothesis, 
first put forward by GARELLI and WATAGHIN. 1.c. (2°), that very energetic 
pions are not absorbed by nuclei. If, on the contrary, energetic pions are 
absorbed in nuclear encounters, it is more natural to extrapolate with a power 
law the generation spectrum of pions, rather than the sea level spectrum of 
muons; if this is done, however, the competition between z-u-decay and 
absorption of pions leads us to expect a steeper decrease of the sea level muon 
intensity with increasing energy. Smaller losses would then account for the 
same intensity-depth curve. 

No direct evidence can be brought forward either in favour or against the 
legitimacy of neglecting nuclear absorption of very energetic pions; however, 
such indirect evidence as exists, (e.g., the inverted temperature effect under- 
ground (3°) (#1), and the anomalous zenithal distribution (°?) (°?), seems rather 
to disprove the hypothesized change of the pions’ behaviour at very high 
energies. 

We are therefore inclined to accept Hayakawa’s opinion and will carry 
out the range evaluation, first, on the basis of «extra-nuclear » losses alone; 
we will then introduce, for the purpose of comparison with future experiments, 
a tentative value for fractional «nuclear » losses, namely 1.0 X10-$ g”!cm?, linear 
interpolation being certainly allowed between the two curves. 


it refers to an inclined direction and has an enormous statistical error; besides, Clay’s 
intensities have been plotted by GEORGE against the depth values given by Cray, 
which were not evaluated on the usual density basis. If the « mass-equivalent » values 
are introduced, the agreement between Clay’s data and the average curve appears 
to improve. Finally recent absolute measurements of intensity by BOLLINGER and 
Cocconr at 1574 hg/em? (private communication) fit rather well George’s curve. 

(2) D. E. Caro, J. K. Parry and H. D. RATHGEBER: Nature, 165, 689 (1950). 

(50) M. Forro: Phys. Rev., 72, 868 (1947). 

(31) J. W. McAnurr: Ph. D. Thesis London (1951 as quoted b. GEORGE, L. c. (16), 
page 420). 

(3?) J. Barnotuy and M. Forro: Zeits. f. Phys., 104, 744 (1937). 

(#3) Cocconr and others, private communication; the effect is quite well established 
by these authors. 
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6. — Chemical composition of the rock. 


Very often in experiments underground the chemical composition of the 
overlaying rock is not perfectly known, while a fairly good estimate of its 
average density is easier to get. 

Collision losses are proportional to Z/A; it is to be expected and easily 
verified that, except for the case of water, no serious error will be introduced 
by putting Z/A — 1/2 (34); but, at high energies, where radiation and pair 
creation losses, which depend on Z?/A, play an important role, the departure 
of single rocks from the average behaviour may be expected to be more serious. 

To investigate this point we first evaluated the average value of Z?/A for 
the average composition of various classes of rocks and we got: a) (Z?/A)ay=5.9 
for igneous rocks; b) (Z?/A)ay = 5.6 for sedimentary rocks (#5). Crystalline 
schist lay obviously in between. 

We then examined some typical rocks in each class and found that the 
single average values of Z?/A ranged from 5.5 (Pegmatite, Elba Island, with 
no iron content) (?¢) to 6.3 (Olivine basalt, Hurtsville, near Sidney, New South 
Wales, with a density of 3.03 g/cm?) (7) for igneous rocks, while departures 
from the average are generally smaller for sedimentary rocks. It is to be 
noted that the above given extreme values are exceptional and refer to rocks 
comparatively rare in the outer crust of the earth. Of course for simple 
rocks departures may occasionally be bigger (e.g. kaolin gives (7?/4),,— 5.0 and 
salt (Z?/A)ay~ 7.0); but simple rocks constitute practically no problem, owing 
to their well known chemical composition. 

It turns out, therefore, that the situation is better than could be expected at 
first thought; as a matter of fact, even if we assume for Z?/A the round value 6 
(which will generally be in excess, especially for sedimentary rocks), it is 
readily seen that the error in the range, up to 500 GeV, will not exceed 2% 
in the worst case. 

For the same typical rocks the average value of Z/A was also evaluated 
and it was found to be always less than 0.5, but with an error which never 
exceeds 1%. 

«Nuclear » losses would be independent of the composition of the rock if 
purely nucleonic absorption were assumed; but, even if we consider these 


(84) At low energies, of the order 10 GeV, of course, the exact knowledge of the 
composition of the rock would be needed for an accurate evaluation of the density 
effect correction, cfr. note (?). 

35) K. Rankama and Tu. G. Samama: Geochemistry (Chicago, 1949), page 222. 
36) Hneiclopedia Italiana, vol. XVII (Roma, 1933), page 724. 
37) B. S. WASHINGTON: Chemical Analyses of igneous rocks (Washington, 1917). 
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losses to be due to interaction with the nucleus as a whole, and assume a 
cross section proportional to the geometrical one, no substantial error will be 
introduced by a purely massic absorption, in view of the other existing uncer- 
tainties on these losses. 

In conclusion we will carry out the range calculation for a standard rock 
(s.r.) defined by A,, = 24, Zay = 12, (Z/A)ay= 0.5, (Z?/A)ay= 6 and we hope 
that the results will apply practically to every rock. Anyway, to take into 
account extreme cases, we will put the results in a form which allows speedy 
correction for any actual value of (Z?/A); besides, the case of water (Z/A=0.555; 
Z?/A = 3.66) will be dealt with separately. 


7. — Range in standard rock and in water. 


On the basis of the above detailed expressions (6), (8), (18) for the va- 
rious types of energy loss, the energy range relation can be evaluated by 
numerical integration. We used the following values of the basic constants: 
ro = 2.82 X10-% cm, me? = 0.511 MeV, N= 6.02 X10%, « = 1/137, _ u-meson 
mass = 210 m (88) and a~ 60. 

The results, for the energy interval 5 to 100 GeV, are collected in Table I, 
with a numerical accuracy which obviously exceeds the actual uncertainty, 
due to the various causes previously discussed. The table also includes the 
« collision range » in s.r., namely the range that would result if collision losses 
were only present: this enables one to evaluate corrections for different values 
of (Z?/A) by simple interpolation. 

Numerical integration would of course do also for higher energies; it is, 
however, expedient, in order to be able to introduce easily any corrections 
that may be required (e.g. to take into account a possibly different value of 
(Z?/A)ay) to put the energy-range relation in analytical form. 

Now it is easily seen that the total energy loss may be written as 


dB 


(19) (D). =A + BBB, 


where A(E) and B(E) are slowly varying functions of #. By assuming them 
to be constant and by putting C — A/B, the integration of (19) gives for the 
range R(E): 


(20) R(E) = R(E,) + a as 


(88) In the above quoted paper (17) we assumed y = 215 m, but this makes no 
significant difference in the numerical coefficients within brackets in formulae (9) and (18). 
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TABLE I. 

| | Range in water hg/cm? Range in standard rock hg/em? 
Energy «extra | «nuclear st «extra « nuclear 
GeV | nuclear» | losses » collision nuclear » losses » 
losses alone | included losses alone losses alone included 

5 | 21.6 21.6 25.0 25.0 25.0 

10 | 42.8 42.8 48.6 48.6 48.6 

| 20 83.6 83.6 93.9 93.9 93.9 

| 30 | 122.8 122.4 138.0 137.0 136.5 

40 161.1 160.3 181.5 178.9 17729 

50 | 198.9 197.4 224.5 220.0 218.2 

60 | 236.1 233.8 267.1 260.3 DENT 

70 272.8 269.7 309.5 300.0 296.3 

80 309.1 305.1 351.6 338.9 334.2 

90 345.1 | 339.9 393.5 377.3 371.4 

100 380.7 374.4 435.3 415.2 407.9 


which has actually been used by Hayakawa and TomonAGAa, l.c. (#), and 
by GEORGE, 1.c. (1%), with suitably chosen values of A and B. 

We have split the energy interval 100-500 GeV in four subintervals 100 GeV 
wide, and determined the relevant values of the constants from formulae (6), 
(8), (18). These values are collected in Table II (energy in GeV, range in 


TABLE 11. 


Auxiliary constants for the evaluation of the range (« extra-nuclear » losses alone). 


E, R(E,) 104 Br, Cz, R(E, + 100) 
standard rock 

100 415 2.680 908.0 768 

200 | 768 2.840 873.1 1 082 

300 | 1 082 2.936 853.9 1 365 

400 1365 3.006 841.1 1 623 

water 

100 | 381 1.635 1 652.3 720 

200 720 1.732 1 588.7 1034 
° 300 | 1034 1.791 1 553.8 1327 

400 1327 1.833 1530.6 1603 
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hg/em?); the subscript Æ, denotes an average value over the 100 GeV sub- 
interval beginning at 2). 

By means of formula (20) an additional «nuclear » loss is also easily taken 
into account; formula (20) could also be tentatively used to extrapolate, 
with suitably chosen values of A and B, the energy range relation beyond 
500 GeV. 

These results, however, would be less reliable on many accounts, as it has 
been made clear in the preceding paragraphs. 

Finally the energy range relation has been plotted in fig. 1 and in fig. 2; 


| GeV. GeV. 
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Fig. 1. Fig. 2. 


upper solid curve refers to s.r., lower solid curve to water, both evaluated 
with «extranuclear » losses alone; dotted curves include instead a fractional 
energy loss of 1.0 x10-* g-1 cm?; interpolation between solid and dotted curves 
will allow to take «nuclear» losses into account, once they will be better 
known. 


In conclusion we wish to thank Prof. E. AmALDI and Prof. S. FRANCHETTI 
for their stimulating interest in the present work and Prof. CocconI for com- 
municating his results prior to publication. 
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RIASSUNTO 


Nessuna delle relazioni fra energia e percorso finora usate per 1 mesoni in roccia 
puö essere considerata soddisfacente quantitativamente. Nel presente lavoro si discu- 
tono le varie difficolta, si migliorano alcune valutazioni approssimate gia adottate da 
altri autori, e infine si valuta il percorso, tenendo conto solo delle perdite per urto, 
radiazione e creazione di coppie. Riteniamo che questa valutazione riduca le incertezze 
esistenti nel confronto fra teoria e esperienza, e permetta di dare sperimentalmente 
dei limiti più stretti alle perdite « nucleari ». A questo scopo la relazione energia-percorso 
è stata valutata anche con un ipotetico valore di 1,0 X 10-6 g-1 cm? per la perdita di 
energia relativa dovuta ad eventi «nucleari ». 


NODE TE CN ECs 


Applicazione del metodo del galleggiante alla misura rapida 
della concentrazione in D,O di soluzioni acquose di acqua pesante. 


E. Cerrar, C. MARCHETTI e M. SILVESTRI 
Laboratori OISE - Milano 


(ricevuto il 21 Aprile 1952) 


Riassunto. — Si descrive un metodo di misura della concentrazione di 
D,O in campioni di acqua, fondato sull’uso del galleggiante. Il metodo 
non richiede taratura con campioni aventi concentrazione nota in D,0. 


1. — Si vuol descrivere l’impiego di un metodo ben noto, quello del galleg- 
giante, nella variante introdotta da SAPIRSTEIN (+), utilizzata perd come metodo 
assoluto, cioè senza l’impiego di campioni di concentrazione nota in D,O, per 
la taratura dell’apparecchio. 

Un campione di acqua da analizzare, avente un volume da 0,5 a 2 cm? 
opportunamente neutralizzato, se necessario, e trattato con qualche me di 
KMnO, viene introdotto nella bolla pit alta di un distillatore a tre stadi, 
costruito in vetro neutro, secondo lo schema di fig. 1. L’introduzione si fa a 
mezzo di siringa in vetro con ago di acciaio inossidabile, dopo di che si applica 
il cono A, ingrassato con grasso da vuoto, e'si fa il vuoto nell’apparato con 
una pompa rotativa munita di trappola a ghiaccio secco. Si fa in modo che 
il degasamento del campione (che dura qualche decina di secondi) avvenga 
senza eccessivo ribollimento, dopo di che l’aspirazione viene esclusa, e l’ap- 
parato, in cui la pressione è di aleuni mm Hg, è pronto per la distillazione. 
Basta chiudere con appositi stringigomma i condotti (1), (3) e (4) ed accen- 
dere la resistenza (1') che riscalda l’acqua, intorno alla prima bolla; il cam- 
pione ne evapora e condensa nella sottostante, che & raffreddata con acqua 
corrente (a circa 15 °C) attraverso il condotto (2). Quando la prima distilla- 
zione è compiuta, si procede con le altre, usando analogamente le resistenza (2') 
e (3’) ed il raffreddamento ad acqua. L’ultima distillazione porta direttamente 
il campione nello stelo (C) in cui sgecciola, attraverso l’ugello (w) avente un 


(4) L. A. SAPIRSTEIN: Journ. Lab. Clin. Med., 35, 793 (1950). 
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diametro di eirca 2 mm, depositandosi nel rigonfiamento del fondo che con- 
tiene un piccolo galleggiante di vetro, fino a raggiungere un livello pari a cirea 
una volta e mezza o due l’altezza del galleggiante stesso. Questo procedimento 
di preparazione del campione non impegna 
pit di mezz’ora o tre quarti d’ora. 


——vuot . 
st La potenza delle resistenze (1'), (2') e 
: (3) puö essere variata da circa 20 a circa 
uscita acqua di = En AUS 
raffreddamento 80 W; esse sono costruite in filo di nichel- 
cromo, immerse direttamente in acqua ed 
ingresso aria secca alimentate a circa 40 V. Si usa la mas- 


sima potenza per scaldare l’ampolla che 
contie le il liquido in evaporazione e la 
minima per tutte quelle a monte di essa. 

Ottenuto il campione puro, si introduce 
aria secca nell’apparato e si distacca lo 


27cm 
MT TE 


ingresso acqua di raffreddamento 


stelo (C) che viene collocato in apposita sede in modo che peschi in un bagno 
termostatico come in fig. 2. Lo stelo pud essere chiuso con cono di vetro quando 
si misurano concentrazioni superiori al 5% in D,O, per evitare inquinamenti 
dovuti all’umidità atmosferica durante il breve tempo di misura. 

I] bagno termostatico & costituito da una vaschetta di vetro a facce all’in- 
circa piane e parallele (che pud essere delle seguenti dimensioni: 22x15 x12 em?) 
quasi completamente riempita di acqua tenuta in forte agitazione dall’agita- 
tore (A), azionato dal motore (M). Un tubo di rame (s), percorso da acqua 
di rubinetto con portata regolabile, consente un raffreddamento dell’acqua 
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della vaschetta, quando necessario, mentre questa pud essere riscaldata me- 
diante Ja resistenza (R) di potenza da 100 a 600 W in nichel-cromo a 40 V 
che vi à immersa direttamente. In tal modo si puö ottenere qualunque tempe- 
ratura compresa fra + 20°C e + 60°C con la desiderata velocita di ascesa 
o di discesa. Accanto allo stelo (C) si trova il termometro (t) che puö essere 
al decimo di grado o al centesimo, del tipo Beckmann. Le temperature souo 
lette a mezzo di cannocchiale, e quindi sul Beckmann é apprezzabile la tem- 
peratura fino al millesimo di grado, il che è veramente sovrabbondante a causa 
di altre fonti di errore che limitano la sensibilità del metodo. Il galleggiante, 
in vetro Pyrex 0 Ignis, è contenuto nel fondo dello stelo che & fatto ad oliva 
in modo che non ci sia azione di parete fra galleggiante e recipiente, e deve 
avere densità apparente uguale a quella dell’acqua naturale ad una tempe- 
ratura compresa all’incirea tra i 25 °C ed i 30 °C. 


2. — La misura si effettua nel modo seguente: quando tutto & disposto 
come in fig. 2, si riscalda, mediante la resistenza (R), l’acqua della vaschetta, 
finchè la densité del campione non sia scesa ad un punto tale, che il galleg- 
giante si immerge fino a toccare il fondo dello stelo; quindi si interrompe il 
riscaldamento e si lascia raffreddare il bagno termostatico per dispersione natu- 
rale o mediante il refrigerante (s) in modo perd, che la velocita di discesa della 
temperatura non sia superiore a 0,2 °Ö/min, perchè la temperatura del cam- 
pione, nello stelo, coincida con i temperatura che si legge sul termometro 
entro la precisione con cui si fa la lettura. 

La densité del campione salira allora gradatamente e, quando essa rag- 
giungerä il valore della densité apparente del galleggiante, quest’ultimo si stac- 
chera dal fondo ed inizierà la salita attraverso il liquido, man mano che la 
densité del mezzo continuerà ad aumentare. L’osservazione del distacco dal 
fondo é fatta o traguardando attraverso la bacinella, a facce pressapoco piane 
e parallele, del bagno, o proiettando un’immagine ingrandita. Si prende nota 
della temperatura t, a cui avviene il distacco (ripetendo, per sicurezza, la iet- 
tura). Si ottiene la Di entro -+ 0,01 °C, e l’attrito di primo distacco 
& praticamente nullo, perché lo stelo si trova in Rane vibrazione dovuta all’agi- 
tazione dell’acqua. Se si ripete la stessa operazione con un nuovo campione si 
otterra la temperatura ¢ di distacco; orbene la differenza t—t, à funzione, 
a parità di ¢,, della concentrazione in D,O del campione. 

Come é noto (?), la relazione che lega la densité di un’acqua alla sua con- 
centrazione æ in D,O è data dall’espressione 2 = aAS/(1— bAs), dove AS = 
= Ô/d; — d,/d; (à è la densità della miscela incognita, d/ la densita dell’acqua 
naturale, d, la densita dell’acqua, H,O, priva di deuterio, il tutto alla stessa 
temperatura t,ed a e b sono opportuni coeffeienti, funzioni della temperatura) 
Se {, € la temperatura di distacco del galleggiante in acqua naturale, la sua 
densita d, sara uguale a quella dell’acqua naturale a quella temperatura (nel- 
Vincertezza di + 0,01 °C), mentre la sua densita d alla temperatura di di- 
stacco ¢ nell’acqua in esame, sara: d = do/(1 + 3/(t—t)), dove 38 é il coeffi- 
ciente di dilatazione cubica del vetro di cui é costruito il galleggiante Ma d 
è pure la densitä della miscela incognita (sempre nei limiti della incertezza 


(?) I. KIRSCHENBAUM: Physical Properties and analysis of heavy water (New York, 
1951), p. 15. 
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corrispondente a + 0,01 °C) e quindi si pud introdurre questa espressione 
nelle formule su scritte, in luogo di 6. Tenendo conto delle espressioni di a 
e di b si ricava: 

: do 1 
(1) a cd;[1 + 3B(t — to) 


= j M, a, ; ode x, 
M, | M, a, ee] 


in eui, riassumendo: 


d, = densita dell’acqua naturale alla temperatura ty, 
di = densità dell’acqua naturale alla temperatura t, 
d, = densitä dell’H,O pura alla temperatura t, 

d, = densità della D,O pura alla temperatura t, 


ce = 0,999984, 
M, = 20,0338, 
M, = 18,0206, 
M, 


— = 1,11172. 
MM 


Data l’omogeneitä della formula, le densitä possono essere sia assolute che 
relative. 

Cosi, per ogni tipo di galleggiante, resta stabilito B, dalla qualita del vetro, 
e da questa e dalla forma la 
sua densità apparente e 
quindi la temperatura 1,, 
percid si possono tabulare i 
valori di w per diverse diffe- 
renze di temperature ¢—t,, 
in modo che, costruendo un 
grafico, siano ricavabili tutti 
i valori intermedi. Assu- 
mendo una nuova tempera- 
tura {, Si ha una nuova curva 
the & ugualmente significa- 
tiva e che serve per galleg- 
gianti che hanno la tempe- 
ratura di distacco in acqua 
naturale in quell’intorno. Se 
perd vogliamo una curva in- 
dipendente dal coefficiente 
di dilatazione del vetro f, 
in modo cioè che sia valida sempre, e che l’influenza di quest’ultimo venga 
introdotta come correzione, si pud tener conto che essendo 3ß(t — t,) sempre 
<1 si ha, con ottima approssimazione d,/(1+ 36(¢—t))= ~ d,{1— 38(t — t)], 
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ed allora la formula (1), con opportune semplificazioni, diviene: 


d d, 

Bary ee 

Nee a re ed; 
= (ie | e& Ww M, LAS 
M, ) M, } cd; M,\ 2) 


In tal modo z risulta funzione di t e t,; per ogni galleggiante (f,) esiste 
una certa funzione f(At) = fı(t,, £) — ffs(t,,t) che fornisce la concentrazione 
in D,0, in funzione della differenza t—t,. La funzione f, dä la concentra- 
zione, supposta nulla la dilatazione del vetro del galleggiante, la seconda, mol- 
tiplicata per f, fornisce il termine correttivo dovuto a questa causa. Per con- 
centrazioni inferiori a circa 0,2% questa correzione puô essere trascurata. 
A scopo indicativo riportiamo in fig. 3 un grafico della funzione f,. Per la fun- 
zione ff, si pud assumere (per galleggianti di pyrex o ignis) l’espressione appros- 
simata (per At< 10°C): ff, % =~ 0,01At. 


SUMMARY 


A float method for determining heavy water concentration in water samples of 
low D,0 content (a few mole percent) is described. This method needs no calibration 
with samples of known D,O content. 
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Rivelatore a scintillazione per particelle «. 


A. MALYICINI 
Laboratori CISE - Milano 


(riceyuto il 21 Aprile 1952) 


Riassunto. — Si descrive un contatore a scintillazione per la rivelazione 
di contaminazioni radioattive di emettitori 4. Sue principali caratteri- 
stiche sono: alta sensibilita, semplicità di costruzione e praticita d’uso. 
Tale contatore & insensibile alle radiazioni X ey e in condizioni di buona 
rivelazione per le particelle « ha un fondo inferiore a 10 impulsi per ora. 


1. — Un contatore a scintillazione per la rivelazione di particelle x consiste 
essenzialmente di uno schermo fluorescente, un tubo fotomoltiplicatore, una 
scatola di protezione dalla luce e un sistema di amplificazione e conteggio (+) (2) (3). 

Abbiamo usato come schermo fluorescente il solfuro di zinco attivato con 
argento, in strati molto sottili; questo fosforo & particolarmente adatto nel 
caso di particelle z le quali hanno un percorso corto. Per strati di spessore 
dell’ordine di una decina di mg em? il solfuro di zinco assorbe completamente 
le particelle pur rimanendo pressochè trasparente alle radiazioni luminose da 
esso emesse. 

Sono stati studiati diversi procedimenti per la preparazione degli strati. 
I migliori risultati sono stati ottenuti facendo sottili pellicole con una miscela 
di ZnS in polvere, formvar e cloroformio nelle seguenti proporzioni: 


Zn D... Der Er > x 2 
FOOT APE EME ee Co. We Ob 
OO ÉOMME RE ML 22 20 


Depositando la soluzione indicata su di un piano levigato e perfettamente 
orizzontale e lasciandola evaporare, si ottengono sottili pellicole fluorescenti. 


() F. H. MazsHaıLı, I. W. GoLpman e A. J. BENNET — A.S.I., 19, 744 (1948). 
(2) J. D. Graves e J. P. Dyson: R.S.I., 20, 560 (1949). 
(?) A. 8. Gotprs, E. R. ROHRER e R. L. Mackin: R.S.I., 21, 554 (1950). 
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Ottimi risultati sono stati pure ottenuti usando polvere microcristallina di 
ZnS semplicemente incollata su cristalplast in modo da ottenere uno spessore 
di ZnS di 10 mg em. 

I «fosfori » deseritti sono insensibili ai raggi X e y. 

Il grande rendimento luminoso del ZnS rende possibile l’impiego di comuni 
fototubi RCA 931 A, nonostante il loro ele- 
vato fondo. A causa del corto percorso delle 
particelle « la protezione del fototubo dalla 
luce in corrispondenza dello strato fluore- 
scente deve essere ottenuta con una parete 
molto sottile, attraversabile dalle particelle «. 

Abbiamo usato a questo scopo dei fogli 
sottili di alluminio in duplice strato. Un solo 
strato di alluminio, a causa della sua poro- 
sita lascia passare la luce in modo da au- 
mentare notevolmente il fondo del fototubo. 

Il circuito di alimentazione del fototubo 
e illustrato in fig. 1. Gli impulsi vengono 
presi dall’ultimo dinodo anziché dall’anodo il 
quale & collegato a terra. Si evita in questo 
modo di rivelare insieme al segnale le scariche 
fra il catodo e l’anodo trai quali esiste una 
differenza di potenziale di circa 1000 V in 


Fig. 1. — Connessioni al fototu- condizioni normali di lavoro. 

bo. Le resistenze sono tutte di L’impulso cosi ottenuto non e sensibil- 

300 kQ. F, fotocatodo; X, elet- mente ridotto in ampiezza ed ha segno in- 

trodo collettore; U, uscita del-  vertito rispetto all’impulso che si ottiene co- 
l’impulso. munemente dall’anodo. 


2. — Sono stati realizzati due tipi di rivelatori. 

Il primo tipo à illustrato nella fig. 2. 

Il fototubo & fissato in una scatola parallelepipeda che porta una finestra 
nella faccia prospiciente al fotocatodo. Su questa finestra viene fissato lo strato 
fluorescente; una guida di luce di alluminio a piramide tronca convoglia i 
fotoni al fotocatodo. 

La protezione dalla luce à realizzata disponendo sopra il fosforo per tutta 
Varea della finestra due strati di alluminio di 5 u di spessore ciascuno. 

L’energia assorbita da questo schermo €, per particelle incidenti perpendi- 
colarmente, di circa 2 MeV. In conseguenza di cid viene ridotto il numero dei 
fotoni emessi dal fosforo; ma tale perdita & compensata dalla riflessione che 
lo strato di Al opera sui fotoni e l’altezza dell’impulso nel fototubo non viene 
sensibilmente ridotta. 

Si è misurata l’efficienza (rapporto fra il numero degli impulsi e il numero 
delle particelle incidenti) per un fascio di particelle « collimate perpendicolar- 
mente allo strato di Ale con energia di qualche MeV, utilizzando una sorgente 
di Po tarata a mezzo di una camera di ionizzazione. Il valore ottenuto è molto 
prossimo a uno. Per & non collimate, l’assorbimento da parte dell’alluminio, 
riduce l’efficienza a circa 1/2. 

L’area sensibile del rivelatore à circa 4x3 em? e si è controllato che l’efhi- 
cienza & pressoche costante su tutta la superficie. 
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L’altezza di un impulso da particelle « di polonio all’uscita del fotomolti- 
plicatore & di 600 mV ed il rapporto segnale/fondo & circa 40. 

Il secondo tipo di rivelatore, é stato ottenuto, incollando, con balsamo del 
Canada, direttamente il solfuro di zinco sul fotomoltiplicatore in corrispon- 
denza del fotocatodo e ricoprendo l’intero fototubo con due fogli di alluminio 
di 5 u di spessore. In questo modo l’area sensibile viene ridotta a 3 x1,6 em2, 


in compenso l’altezza dell’impulso 
viene piu che raddoppiata in con- 
fronto a quella ottenuta con il tipo 
1°, per la miglior raccolta dei fotoni 
sul fotocatodo. Si hanno impulsi alti 
1,3 V all’uscita del fotomoltiplica- 
tore con una tensione d’alimentazione 
di 1000 V. Il rapporto segnale/fondo 
è circa 80. 

Tale rapporto varia notevolmente 
con il potenziale del fotocatodo, ed 
aumenta quando tale potenziale di- 
minuisce. Facendo lavorare il foto- 
tubo a 650 V con il 2° tipo di conta- 
tore si ha un valore di questo rap- 
porto di circa 120. La variazione del 
rapporto segnale/fondo in funzione 


IRRE 


agree 


à 
a 
= 


Foto- 

moltiplicatore (931)A. 2. Guida di luce di 

3. Schermo di ZnS. 4. Fine- 

stra sottile di alluminio. 5. Reticella me- 
tallica protettiva. 


Fig. 2. — Sezione del rivelatore. 1. 


alluminio. 


della tensione dipende anche dalla corrente di fondo dovuta all’effetto di ten- 
sione del vetro del fotomoltiplicatore (4). 


Fig. 3. — Rivelatore di particelle «. 


(4) M. G. TayLor: Journ. de Phys. et le Rad., s. VIII, 10, 255 (1949). 
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I due contatori schermati nel modo descritto sono insensibili alla luce visi- 
bile, come & stato verificato esponendoli alla luce di una lampada di 100 W 
accesa à pochi cm di distanza. 

In entrambi i tipi di rivelatori abbiamo usato una RC di ingresso di 10? 5; 
si ottiene cosi un potere risolutivo sufficiente allo scopo. 

Data la notevole ampiezza degli impulsi, per la loro registrazione è suffi- 
ciente un amplificatore a guadagno 50-100. Usando un amplificatore del genere 
e un discriminatore con un taglio corrispondente a condizioni di buon funzio- 
namento per il conteggio di particelle x di qualche MeV, il fondo & inferiore 
a 10 impulsi per ora. In queste condizioni, variazioni della tensione d’alimen- 
tazione dell’ordine del 5% non influenzano sensibilmente il conteggio. 


3. — L’apparecchio realizzato in questi laboratori, nella forma indicata 
dalla fotografia fig. 3, e seguendo lo schema di fig. 2, s’e dimostrato efficace 
monitore in manipolazioni intorno al polonio, il quale, come & noto, dimostra 
un elevato grado di tossicita. 


SUMMARY 


The A. describes scintillation counter for the detection of radioactive contamina- 
tions of « particle emitters. Its main features are: high sensibility, simple building and 
easy use. This counter is insensible to X and y rays; under conditions of good « particle 
detection, its background is less than 10 pulses per hour. 


LETTERE ALLA REDAZIONE 
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Hydrogen Isotopes in Steam Wells. 


G. Boato e G. CARERI 


Istituto di Fisica dell’ Universita, 


'entro di Studio per la Fisica Nucleare del C.N.R. - Roma 


G. G. VOLPI 


Istituto di Chimica dell’ Universita - Roma 


(ricevuto il 15 Aprile 1952) 


As far as we know the deuterium 
content in natural gases seems not to have 
been much considered (1). Following a 
suggestion of the Larderello Corp., we 
have made a rough analysis of the hydro- 
gen isotopes of the steam from the Sof- 
fioni (?), and of the elementary hydrogen 
which is also contained in these steam 
wells, in a quantity about 1/1 000 or less 
of that of the steam. At present there are 
no clear ideas about the origin of this 
hydrogen, which is always present in 
natural and volcanic gases. 

To measure the deuterium content in 
the steam, hydrogen has been prepared 
from it by complete reduction of water 
samples with metallic zine at 350 °C. 
Some samples artificially enriched by 


@) A. F. Scott: Science, 79, 565 (1934), 
and L. SHAMOYSKY and N. KAPUSTINSKAYA: 
Acta Physicochim. U.R.S.S., 7, 797 (1937) ana- 
lyzed some water samples collected in the inside 
of the earth at a depth of about 1500 meters, 
finding no significant variation in the heavy 
water content compared with that of tap water. 

(2) For a brief description of these «Sof- 
fioni» of Larderello and bibliographic details 
about them, see our previous note on the argon 
isotopes in: Nuovo Cim-nto, 9, 44 (1952). 


1/10000 part of deuterium have been 
also prepared to test the whole method 
of analysis. A Nier-type mass spectro- 
meter has been used, operated with ma- 
gnetic scanning. The HD/H, ratio has 
been evaluated in the usual way, from 
the extrapolated (HD+— HF) ratio at 
zero H, pressure. The results are as 
follows. The analysis of the water from 
the pipe collecting the condensed steam 
of all the wells, showed an enrichment 
factor (%) 1 + 0.1; the same value has 
been found also fer the steam of a single 
well. The analysis of the hydrogen from 
different wells, showed enrichment fa- 
ctors so near among themselves, to allow 
us to considere a common value for all 
the sources. This value, averaged over 
five steam wells from different places, 
was 0.79 + 0.05. One sees that the 
deuterium content of the steam is not 
different from that of ordinary tap wa- 


‘ ter; but on the other hand the deuterium 


content of the elementary hydrogen is 
certainly rather less. A possible expla- 


[DIA] in the sample 
[D]/[H] in tap water ~ 


(*) Defined as 
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nation of this variation could be found 
in a diffusion separation of the gas, due 
to some porous material inside the earth; 
as a matter of fact a single diffusion se- 
paration could give an enrichment factor 
of (2/3)42 = 0.82. 

An alternative explanation seems also 
to be possible if one thinks of the ex- 
change reaction H,0 + HD = HDO + H,, 
which can be catalyzed inside the earth, 
before the steam can reach the atmo- 
sphere. 

Since the equilibrium constant of 
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this reaction is >1, the deuterium 
must concentrate in the water and leave 
the elementary hydrogen poorer in it. 
Of course one does not observe the en- 
richment of the steam, due to the very 
small amount of hydrogen which is con- 
tained in it. Furthermore if one attempts 
to calculate the temperature at which the 
equilibrium is established from the ob- 
served enrichment factor, one gets a tem- 
perature ranging around 600 + 100 °C, 
which is quite a reasonable one and con- 
sistent with the geological data. 
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Sulle frequenze di vibrazione delle molecole biatomiche 
| formate con atomi del 4°, 5°, 6°, 7° gruppo. 
. 
G. LOVERA 
| Istituto di Fisica dell’ Universita - Modena 
| (ricevuto il 14 mageio 1952) 
| 
| 


La formula, dedotta in una nota precedente (1) 


~ SITE M. : 
(1) » = 1979,7 |/ == 6m-1 (r, 0 AS u In uma): 
ru 


che fornisce i numeri d'onde » delle vibrazioni molecolari, & stata applicata alle 


. molecole biatomiche formate con atomi del 4°, 5°, 6°, 7° gruppo del sistema perio- 
dico. I valori di A’, gid riportati in parte nella suddetta nota, sono i seguenti: 


2 2,5 3 3 
5 Ki 2,146, 2,594, 2,985, 3.345, 3,670, 3,977, 
à nun = 4 4,5 D DD 6 
| RI 4,255 5.431, 4,781, 5.033, 5,270, . 


| A differenza dal caso, trattato nella precedente nota, dei metalli alcalini, per i 
| quali si & posto n = 1, nei caleoli presenti, invece, per le molecole biatomiche i cui 
| due atomi appartengono, nel sistema periodico, allo stesso periodo, si & preso per n 
il numero d’ordine del periodo stesso; nei casi, invece, in cui i due atomi fanno parte 
di periodi diversi, si é posto n uguale alla media dei numeri d’ordine di questi due 
periodi. A v si é dato il valore 2 per le molecole formate con atomi del 4°, 59, 6° gruppo, 
Ê ed il valore 1,6 per le molecole aventi uno od entrambi gli atomi appartenenti al 
7° gruppo. 
Nelle tabelle I e II sono posti a confronto i numeri d'onde » sperimentali con 
quelli calcolati mediante la (1) nelle eondizioni ora precisate. I dati sperimentah 
di r, e di » sono tratti dalle tabelle di HERZBERG (2), tranne i valori contrassegnati 
con (G), che sono dovuti a GUGGENHEIMER (3). 


oN 


() G. Lovers: Nuovo Cimento, 9, 442 (1952). 
3 -(2) G. HERZBERG: Molecular spectra and molecular structure — I: Spectra of diatomic molecules 
+ (New York. 1950). 
(3) K. M. GUGGENHEIMER: Proc. Phys. Soc, 58, 456 (1946). 
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A parte qualche scarto eccezionalmente elevato, l’accordo tra i dati sperimentali 
e quelli calcolati & soddisfacente, anche per le molecole formate con atomi del 4° 
e 6° gruppo, per le quali era gid stata proposta in precedenza un’altra formula, più 
accurata. L’errore medio percentuale, media dei valori assoluti degli errori percen- 
tuali singoli, risulta del 4,5%. 


TABELLA I. 


Numeri d’onde di vibrazione delle molecole formate con atomi del 4°, 5°, 6° gruppo. 


| Mole- 1, y (cm) Diff. ‚Mole- 1 re I (em!) Diff. 
cola | ” (A) 088. cale. % | cola x (A) oss | cal. | % 
| = Fire | 
GRUPPO 49 | | 
LG, 2 | 1,3117 |1641,3 | 1858,8 | +27,6 | | 
GRUPPO 5° | | 
N, 2 |1,094 |2359,6 |2259,1| —4,3 || Bi, 6 |2,92(G) | 172,7 1782| +3,2 
Br 3 1,894 | 780,4| 739,2 | —5,3 || NP | 2,5 | 1,491 1337,2 | 1280,7 | — 4,2 
As, 4 |2,19(G) | 429,4 411,8 | — 4,1 | AsN | 3 |1,65(G) | 1068,0 | 1041,9 | — 2,4 
Sb, 5 | 2,57(G) | 269,8| 269,3 | — 0,2 || SbBi| 5,5 | 2,73(G) | 220,0| 224,56) 2,5 
a SE RE | | | 
—e “ | _ —. 
| Grurpro 6° | | 
| | 
0, 2 |1,2074 |1580,4 | 1832,0 | +15,4 | SO | 2,5 | 1,4933 | 1123,7 | 1216,3 | +8,1 
| Se 3 | 1,889 725,7 | 730,7 | +0,7 | SeO | 3 |1,61(G) | 907,1 | 1018,0 | +12,2 
| Se, A 2,157 391,8 407,8] +4,1 | TeO | 3,5 | 1,82(G) | 796,0 852,9) +7,1 
| Te, 5 2,59 251 259,7 | +3,5 
| GRUPPI 40 e 50 | 
CN 2 1,1718 |2068,7 2121,2 | 19,5 E.SIN a E2,5121°5718111151.700.12028: PART 
CB 12,51] 1,562) 119397) 1261,01 42077 
| re a3 uid Be = La fe el = =| 
| GRUPPI 49 e 60 | 
| 
| co 2 | 1,128 2168,2 | 2180,4! + 0,6 || GeTe| 4,5 | 2,39(G) | 323,41 332,3 | + 2,8 
CS 2,5 | 1,536 1285,1 | 1287,7| +0,2 | SnO | 3,5 | 1,838 822.5 | 843,59) 22,6 
CSe | 3 | 1,735(G)| 1036,0 | 1028,2 | — 0,8 | SnS | 4 | 2,27 488,3 | 475,2 | — 2,7 
S10) 1-255) 1,010. =|71949,0.1 12235221) 21,5 | SnSe| 4,5 | 2,37(G) | 331,2) 335,4| +1,3 
Sis 3 | 1,929 749,5 | 732,9| —2,2 | SnTe| 5 | 2,57(G) | 259,5 267,9) + 3,2 
SiSe | 3,5 | 2,14(G) | 580,0! 554,9 — 4,3 | PbO | 4 | 1,922 721,8 | 786,3| +8,9 
SiTe 4 |2,34(G) | 481,2 476,3 | — 1,0 | PbS | 4,5 | 2,395 428,1 | 431,9| +0,9 
| GeO | 3 | 1,651 985,7 | 958,8| 0,0 | PbSel 5 | 2,53(G) | 277,4 | 287,1 | + 3,5 
| GeS | 3,5 | 2,06(G) | 575,8 | 565,2 | — 1,8 | PbTe| 5,5 | 2,73(G) | 211,8 | 221,6 | +4,6 
| GeSe| 4 |2,19(G) | 406,8 | 406,6 | — 0,0, | | 
| GRUPPI 5° e 6° | 
| NO | 2 | 1,1508 |1904,0 | 2027,6| +6,5 | PO | 2,5 | 1,4488 |1230,6 | 1278,5 | +3,5 
| NS | 2,5 1,50(G) | 1220,0 | 1263,1 | +3,5 | AsO | 3 |1,61(G) | 967,4 |1022,6 | +5,7 


SULLE FREQUENZE DI VIBRAZIONE DELLE MOLECOLE BIATOMICHE ECC. 


Numeri d'onde di vibrazione delle molecole contenenti atomi 


TABELLA II. 


Mole- | me 
cola | d (À) 
nier 


GRUPPO 7° 


Boa, 2") 13485 
Cl, 3 | 1,998 
Br, | 4 | 2,2836 
I |). | 2,6666 
CIF | 2,5 | 1,6281 


1032,8 | +15,7 


GRUPPI 7° e 4° 


COL | 2,5 | 1,73(G) 
SiF | 2,5 N 1,603 | 
Sic] | 3 | 2,00(G) | 
| SiBr | 3,5 | 2,15(G) | 
| GeCl | 3,5 | 2,08(G) | 
| GeBr| 4 | 2,29(G) 


| 
| 


GRUPPI 7° e 5° 


| NBr | 3 |1,91(G) | 
| SbF | 3,5 | 1,82(G) | 


| 
| 


| BiCl | 4,5 | 2,48(G) 


Un accordo migliore con i dati sperimentali si potrebbe ottenere scegliendo, per 
ogni gruppo di molecole, un opportuno valore di v, il che non appare ingiustifi- 


. [Mote- | 
cola | 2 
BrF 3 
BrCl | 3,5 
IC] 4 
[Br | 4,5 


SnE | 3,5 


SnCl| 4 


SnBr 4,5 


PbF | 4 


PbC1 104,5 


PbBr 5 


| 2,60(( 


del 7° gruppo. 


(A) 


2,43(6) 


| 2,63(G) 


2,796) 


>? (em!) 


088. | cale. 


671 665,8 
430 407,1 
384,2 353,0 
268,4 247,5 


582,9 | 635,6 
352,5 | 355,7 
247,7 | 255,9 | 
507,2 | 549,6 
303,8 | 325,3 
207,5 | 218,5 


209,3 | 214,5 
163,9 171,6 


cato, dato il significato fisico attribuito a v nella deduzione della formula (1). 


IL NUOVO CIMENTO Vou. IX, N. 6 1 Giugno 1952 


On a question of nuelear spectroscopy. 


A. GAMBA and M. VERDE 


Istituto di Fisica dell Universita - Torino 


(ricevuto il 16 Maggio 1952) 


After the success of the central model in the study of nuclei, some questions 
have arisen as to the numbering of the states corresponding to given configurations. 
One may ask, for instance, which are the multiplicities of the allowed values of 
the spin and isotopie spin corresponding to a given symmetry class (Wigner’s 


nh 


super- 
multiplet) for an n-nucleon system. This problem has been solved by a general 
formula of Racan (*). 

Another question comes in, when considering the multiplicities N, of the angular 
momenta J associated to a given term of » equivalent nucleons. No general for- 
mula has been given so far for this problem. In atomic spectroscopy both the group 
theoretical and the SLATER (?) method are used. They, however, require that each 
case be considered apart. Recently extensive tables have been published by 
Jaun 6) and FLowers (4) for p, d, f equivalent nucleons, and, in the case of j-j 
coupling, for j = 3/2, 5/2, 7/2, by using a recurrent procedure developed by Janx (°). 

A complete solution of the problem for j=1 may already be found in a for- 
mula due to Racan (ref, (1), formula (20)), who, hewever, has not given the physical 
interpretation of his relation. The method of Racan, which actually consists in 
solving a difference equation, becomes quite impractical for j > 1. 

In this letter we give some general formulae, which solve the problem for 
n == 3,4 and arbitrary j. These formulae are of interest because one has often 
to deal with nuclei with few equivalent nucleons of high angular momenta cutside 
a closed shell. The formulse are valid both in the case of orbital momenta 
(()" states, LS eowpling) and in the case, more significant for nuclei, of total 
angular momenta ((j)® states, j-j coupling). 

The procedure is essentially the group theoretical one, as outlined for the atomic 
cas» by Wer (5). We must write down the irreducible characteristic X of the 


unitary group associated to the partition {ff}, where fi, fy, fg... are the 


(NG, Racamı Rev. Mod. Phys. 21, 494 (1949). 

@) J. ©. Suarar: Phys. Rev, 38, 1109 (1951). 

(9 H. A. Jaan: Proc. Roy. Soc., A, 201, 516 (1950): 205, 192 (1951). 

() B. A. FLowers: Proc. Roy. Soc. A, 210, 497 (1952); 212, 248 (1952). 

(6) H. Wry: The theory of groups and quantum mechanirs (Methuen, 1931), p. 374. 
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numbers of squares in each row of the corresponding Young pattern, One has: 


a ye Ai sfatsn.(Urtaly +) itt eh EN e23\ ta an — EN 
ö 14 Qt: 3... | | es Breas: 


alt... Re es — 1 


where %,,1,.(totiis ++) 18 the character of the irreducible representation of the sym- 
metric permutation group associated with the same Young pattern for the class 
formed by all the permutations with 7, first-order cycles, à, second order cycles, etc., 
and the sum is meant over all such classes. 

The numbers N, are then the coefficients in the development: 


EH GT 


X(e) = >, 


e— 1 
For the particular cases n = 3, 4 the following formulae are easily established: 


RES; A= 3j —J 


1+(—1)4 A +3 A +2 4 +1 
(3){4+1}+ 37(11) ae + 22000 {| 2 |+| = 


— 
N 


for 37 >J >}, 
si | Ze) 
71a {27 +1} + 3711) 1) x a 3 a 
3 Al Hey asi 
+ 27(001) {|* : |-2| u |+ | = |} for j>J>0, 


TELE B = 4j —J ; C = 2 —J 


B? + 3B +2 
4) — 


a 


‘iB 3] Ree B+4 
+ | - | | = oz | 2 


2 2 


B+2 Bene 
+ 6x(21) | | + 3x(02)(— 1)2 "| x | di 


4 


PP 


for 4j >J > 2j, 


= 8jJ + 4j SIEHT 2 C 
Sig ee ent (— 1) a+} + 
B 2 B 3 BE 
+ 3700 | en [+s sion | 3 = =|- 5 2 = 2 
O+2) fo+1 i B+4) .[B+3] [B+2 
| =) con {JP | | ; | | = |} 
| for 2j>4J >0, 
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where [x] means the greatest integer contained in x for 2 >0, and [x] —0 for » <0. 
The extension of these formulae to cases with n > 4 is straightforward in prin- 
ciple, although the calculations become rather tedious. We have not succeeded in 
giving a formula valid for arbitrary n. 
For the reader’s convenience we give the tables of characters of the symmetric 
permutation group for n = 3,4. For n> 4 extensive tables may be found in the 


literature (8), 


N Chefs} | | Ets | | | | 
. u 15) dl an) | Ge 3 {4} | {31} | {22} {211} (1111) 
(ees SS | | | 161029302 | | | | 
> | Be | ie t= aa 
| Pr) | 
2 ee) 1 4 1 3 2 3 1 
e | 2 ee 1 Ge 
LL TE A | 02 LT DEA 1 
I IE ee — ans 
: | | 101 i Owl at 0 1 
| | | ip tr att beeen 


(8) F. D. MURNAGHAN: The theory of group representation (Baltimore, 1938). 
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LIBRI RICEVUTI E RECENSIONI 


LÖSCH-SCHOBLIK — Die Fakultät und 
verwandte Funktionen. Teubner, 
Leipzig, 1951, pagg. 205. 


Quest’opera, come dichiara l’Autore, 
F. Lösch, che l’ha compiuta sulla base 
di un manoscritto lasciato da F. ScHo- 
BLIK, « vuol essere un’introduzione alla 
teoria della funzione fattoriale (funzione 
Gamma) e connesse, la quale accanto 
agli interessi del matematico tenga conto 
anche di quelli del fisico e dell’ingegnere ». 

Questa doppia finalità € tenuta effet- 
tivamente sempre presente, anche nei 
capitoli essenzialmente teorici, cosieche 
di ogni nuova espressione viene posto in 
rilievo l’aspetto applicativo. Che cid non 
vada a detrimento della compiutezza e 
del rigore scientifico, e lasci nello stesso 
tempo del tutto unitario e corente lo 
svolgimento, & certo uno dei prineipali 
pregi del lavoro. 

Non & inutile avvertire che questo 
libro non vuol sostituire aleun manuale 
di caleolo numerico: per quanto riguarda 
quest’ultimo, invero, l’A. rimanda alle 
note tavole di Jahnke-Emde, unifor- 
mando a tale scopo le proprie notazioni 
con quelle ivi usate. 

Nel primo capitolo, dei tre che costi- 
tuiscono l’Opera, vengono considerate le 
varie rappresentazioni della funzione fat- 
toriale J’ (come prodotto infinito, come 
serie, come integrale) nonché questioni 


di calcolo differenziale, integrale e fun- 
zionale che con la J’ si connettono. 

Nel secondo viene considerata la co- 
siddetta funzione fattoriale incompleta. 
che si pud ottenere dall’integrale di Eu- 


co 


lero di seconda specie 2! = | et; dt con 


0 
la decomposizione di Schlömilch 


£ 2 
21 — | edt + jet dt. 
0 e 

Gli integrali a secondo membro ven- 
gono studiati sia come funzioni di g (per 
un dato +) che come funzioni di 2 (per 
un dato g). In particolare il primo dei 
due modi porta a considerare, per deter- 
minati valori di +, importanti funzioni 
come le integral-esponenziali (integral- 
seno e integralcoseno), la funzione degli 
errori, e infine i polinomi di Hermite. 

Il terzo capitolo & dedicato alle appli- 
cazioni: l'A. tocca campi svariatissimi, 
dalla propagazione del calore alla siste- 
mazione dei posti a sedere nei teatri, dai 
vari tipi di radiazioni alle rendite. 

Questo cenno agli argomenti dell’ul- 
tima parte del libro ci induce a segnalare. 
anche per le precedenti, la riechezza di 
contenuto, che la concisione dell’A. è 
riuscita a condensare in modesta mole. 


Dario FÜRST 


ej 


Direttore responsabile: G. POLVANI _ Tipografia Compositori 
Ha 
torehi il 23-V-1952 


Questo fascicolo & stato licen: 


